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Abstract 

In this paper we consider parafermionic Liouville field theory. We study integral represen- 
tations of three-point correlation functions and develop a method allowing us to compute them 
exactly. In particular, we evaluate the generalization of Selberg integral obtained by insertion of 
parafermionic polynomial. Our result is justified by different approach based on dual representation 
of parafermionic Liouville field theory described by three-exponential model. 
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1 Introduction 



Non-rational conformal field theories attract a lot of attention in contemporary theoretical physics since 
Polyakov's suggestion to study strings in non-critical dimension pp. In Polyakov's approach metric on 
a surface of the string due to conformal anomaly become dynamical and its behavior is described by 
Liouville field theory which is the most well known and widely studied in the literature example of 
non-rational CFT [2H5] . An important question which one usually asks about any CFT is the complete 
set of correlation functions of local fields. The conformal symmetry of the theory plus the hypothesis 
about the operator algebra help to answer this question very much In particular, the problem can 
be reduced to the problem of studying of the operator algebra of so called primary fields. In [6J the 
set of conformal theories so called minimal models has been constructed. The main property of these 
theories is a finite number of primary fields. Contrary, the main "disadvantage" of non-rational CFT 
- infinite (even uncountable) number of primary fields was probably the main reason why the progress 
in this direction has been postponed. In [2HHE] the three-point correlation functions or equivalently 
the structure constants of the operator algebra of primary fields where found exactly (DOZZ formula). 
Conceptually, this important result can be viewed as a complete solution of Liouville field theory because 
it allows to find (in principle) any multipoint correlation function. Another piece of information about 
the theory which is needed for the construction of all correlation functions is encoded in the conformal 
blocks which are objects defined only by the conformal symmetry. We note that the conformal blocks 
with general external and intermediate conformal dimensions appear essentially only in non-rational 
CFT. Despite the fact that there is a straightforward algebraic procedure allowing to construct the 
conformal blocks their explicit determination is a tedious problem which become actual again since the 
discovery of the connection between non-rational CFT and four-dimensional M = 2 supersymmetric 
gauge theories in recent paper of Alday, Gaiotto and Tachikawa [S] (AGT relation). Probably the best 
way to compute the expansion of the conformal blocks is the recursion formula suggested by Alyosha 
Zamolodchikov [9j (we note that AGT relation also gives very effective algorithm allowing to compute 
the conformal blocks). 

So far, only Liouville field theory as well as its sl(n) generalizations called Toda field theories [T0TH2] 
were studied in the context of AGT relation. However, we still have many examples of non-rational 
CFT's whose role in gauge theory-non-rational CFT duality is unclear at present. One possible candi- 
date is the supersymmetric generalization of Liouville field theory which appears essentially in Polyakov's 
approach applied to superstring in non-critical dimension [13] . This theory is now studied in the same 
extent as ordinary bosonic LFT. In particular, the supersymmetric generalization of DOZZ formula is 
known [T4"1 [T5] as well as the recursion formula for the conformal blocks [T6HT8] . Possible generalization 
of superstring theory called fractional superstring theory was suggested in [B5] which is roughly speaking 
based on idea to substitute ordinary fermions by Z N parafermions [2D] . 

Theories with parafermionic symmetry were studied in different contexts (see papers citing [20J) and 
their massive deformations in [2T1425] . We give an overview of Zn parafermionic CFT in appendix lAl 
Probably, the most promising application of this CFT appears in studying of non-Abelian Fractional 
Quantum Hall Effect So called Read-Rezayi state constructed in gives trial ground state 
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wavefunction for the system of charged particles on a sphere of certain radius in a uniform radial magnetic 
field at Landau level filling N/ (2 + NM) where M is an integer number. Up to some standard factors 
the (N, 2) Read-Rezayi state is proportional to the multipoint correlation function of parafermionic 
currents whose non-trivial part is given by a symmetric homogeneous polynomial vanishing when any 
N + 1 particles collide. We refer interested reader to recent papers [271431] for more details and list 
of all relevant references. In this paper we study these and even more general polynomials which can 
be identified due to [21] with three quasihole excitations over the quantum Hall ground state in great 
details. In particular, we give an explicit representation for them in terms of symmetrization of certain 
basic "monomials" and compute exactly generalized Selberg integral involving these polynomials (see 
section [3]). We also compute the generalization of Mehta integral involving parafermionic polynomial. 

Fractional superstring theory was investigated in number of papers [32-36J and recently in [371138] 
in the context of matrix models. However, many important questions about fractional superstings are 
still unanswered. One of them is the formulation of non-critical fractional superstring theory within 
Polyakov's approach similar to [UH3]. This question requires deeper understanding of the nature of 
parafermionic symmetry, but it is clear that Zn parafermionic analog of Liouville field theory plays 
here an important role. Parafermionic Liouville field theory is defined by two-dimensional Lagrangian 
density [39J 

£ = £ PF + ^ (d a if) 2 - //iVMe 26 *, (1.1) 

where we denote schematically the Lagrangian of Zn parafermionic CFT by £pp [20] (see also appendix 
Rjl . ip is bosonic field and \I/ is the first parafermionic current. Parameter \i in f ll.ip is called cosmo logical 
constant, while b is coupling constant. Last term in f ll.ip can be viewed as the simplest possible marginal 
deformation of two theories: free bosonic theory and Z^ parafermionic CFT. For N — 1 the resulting 
theory coincides with ordinary Liouville field theory [HE] and for N = 2 with supersymmetric Liouville 
field theory [14HT6] . For arbitrary values of the parameter fi the Lagrangian ( 11.11) defines CFT with 
central charge 

c=— 2 +^, where Q = b+-. (1.2) 

The symmetry algebra of the theory is actually wider and generated by two holomorphic currents of spin 
two (stress energy tensor) and non-local current of spin for N > 1. For N = 2 this algebra coincides 
with Neveu-Shwarz algebra while for N > 2 it has more complicated structure^. For our purposes we 
found inconvenient to use standard "bootstrap" machinery based on using of Ward identities an so on [6] 
and prefer to consider the theory fll.ip as a deformation of the tensor product of two theories free bosonic 
field theory and Zj^ parafermionic CFT. It means that we will consider all correlation functions in our 
model only if the total charge takes special values and they can be represented by finite-dimensional 
integral of the product of correlation functions in free bosonic theory and parafermionic CFT (see 
section [2] for details). The expression for general values of the total charge can be obtained by means 
of "analytical" continuation (see section [2]). This approach is not mathematically rigorous and the 
obtained result has to be checked by independent calculations. In section H] we consider the theory fll.ip 
using the dual representation by the three exponential model which is introduced there and develop a 
method allowing to study correlation functions without any restrictions on the total charge. 
The primary fields in the theory (11.11) are given by the exponential fields 

V™ = a k e 2a \ (1.3) 
1 For N — 4 this algebra can be realized as spin 4/3 parafermionic algebra [40] (see [34TI36J ) . 
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where cx^ are the primary fields in parafermionic CFT also called order fields (see appendix EJ. The 
fields ( II. 3p are the spinless fields with conformal dimensions 



AW(a)=AW(g)= k{N - k) a(Q - a) 



2N(N + 2) N ' 

In this paper we consider in details the three-point correlation functions of the primary fields (II .3p in 
parafermionic LFT and obtain generalization of DOZZ formula [2]-[I] in this case. The paper organized 
as follows. In section [2] we formulate the problem and propose an answer for the three-point correlation 
function of exponential fields in our model. We also describe general method used in this paper which 
leads to the problem of computation of some multi-dimensional integrals analogous to those considered 
in jH] but with insertion of certain polynomials (we call them parafermionic polynomials). In section 
[3] (which can be considered as a main novel part of our paper) we compute these integrals using the 
method adopted from the seminal paper by Selberg [12] • In section H] we consider dual representation of 
parafermionic LFT given by three exponential model and obtain an analytic expression for the three- 
point correlation function which agrees with the result obtained in section [31 Section H] can be considered 
as an independent check of our results. In this part of the paper we use the notations historically related 
with the three-exponential model and they intersect with the notations used in the rest of the text (the 
same is for appendix [D]) . We hope that this fact will not confuse the reader. In section Owe consider 
some important questions which were not discussed in the main body of the paper. In appendix [X] 
we give a review of parafermionic CFT [20] . In appendices [B] and O we give proofs of two important 
propositions essentially used in section [3j In appendix [D] we consider minisuperspace approximation for 
the three-exponential CFT and discuss its sigma-model representation. 

2 Three-point correlation function in parafermionic LFT 

Our goal is to find the three-point correlation function of the fields (II. 3p 

C klMM (a l7 a 2 ,a 3 ) d = f <y<£Ho)V£\l)V<*\<x>)). (2.1) 

One can show (see below) that due to the symmetry of parafermionic CFT the three-point correlation 
function (12. ip is non-zero only if one of the following conditions is satisfied^] 

k 1 + k 2 + k 3 = 2k, < kj < k < N, (2.2a) 

or 

k 1 + k 2 + k 3 = N + 2k, < k < kj < N. (2.2b) 

In this paper we give an analytic expression for the three-point correlation function (I2.ip . In order to 
present an answer it is convenient to define functions 



C kM ( ai ,a 2 ,a 3 ) 



71 m \~n) n 



Q-a 
b 



p(ki,k 2 ,k 3 )x 
Tl N \b)T[ N i \2a 1 )T^\2a 2 )Ti N 3 \2a 3 ) 



(2.3a) 

We note that for even N both conditions (|2.2j) can be satisfied at the same time. We discuss this phenomenon below. 



T fc Ar) («i + "2 + ol 3 - Q)T^} ki (a 2 + a 3 - a^T^l^ + a 3 - a 2 )TfJ k3 (ai + a 2 - a 3 ) 



and 

C fcl ' fc2 ' fc3 ( ai ,a 2 ,a 3 ) 



\ ,_26« 

717/7 ( — J 6 » 



Q-a 



p{h, h, h)x 
Tr ) (&)Tir ) (2« 1 )Tlf ) (2« 2 )T^(2a 3 ) 



Tfc (ai + a 2 + a 3 - Q)T^J fc _ fa (a 2 + a 3 - ai)T]JJ fc _ fa (ai + a 3 - a 2 )T ( ff} k _ k3 {a 1 + a 2 - a 3 ) 



(2.3b) 



where a = a\ + a 2 + a 3 and integer parameters fcj and k are supposed to be related by (I2.2a[) for the 
function (j 2. 3 aft and by (I2.2bl) for the function (12.3bj) . The factors p(k±, k 2 , k 3 ) and p(k\, k 2 , k 3 ) in f)2.3ap 
and (12.3bj) are functions independent on parameters a and b. They are explicitly given by 



p(ki,k 2 ,k 3 ) 

p(h, k 2 , k 3 ) 
where 



7 



7 



iV + 2 
1 

N + 2 



n-, 



iV - kj + 1 
iV + 2 

N-kj + 1 
N + 2 



H(k) = l[ 1 



#(iV - k)H(k - ki)H(k - k 2 )H{k - k 3 
H(k 1 )H(k 2 )H(k 3 ) 

H(k)H(h - k)H(k 2 - k)H(k 3 - k) 
H(h)H(k 2 )H(k 3 ) ' 



J 



(2.4) 



N + 2 

Throughout this paper we use the following standard notation 

T(x) 



7(x) 



r(i -x 

In eqs (I2.3P function T k N \x) = T k N \x,b) is defined as^| 

'x + kb- 1 + (j - l)Q 



(2.5) 



JV-* 



Tf)(a;)=n T 



i=i 



N 

n t 

j=N-k+l 



x + (k- N)b~ l + (j - 1)Q 



(2.6) 



N / 11 V N 

/ j=N-k+l v 

where T(x) is usual T-function defined in [I]. It is convenient to supply definition (12 .6ft with 

which is self-consistent for k = as one can readily see from (12.61) . Function T k N \x) obeys several 
important properties (which follow of course from analogous properties [I] for T function). 



Reflection: 



T k N) (Q-x) = T 



(AT) 
N-k 



. X 



(2.7a) 



^Usually the b dependence of the function t£ (a:, b) is hidden for shortness. 
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Shift: 



r(")(x + 6 ) = 6-^ 7 f^±^)rK(x 



N 



A 7 



Duality: 



Tf ) (x,6) = T^ fc (x,r 1 ). (2.7c) 



We note that due to the relation ( I2.7bl) all function Tjjj. (x) can be obtained from the function T^ A ,!| 
given by integral representation 



x 



log T-V) = f f (f - x)* «- - ^^t^-f ) f ■ (2-8) 

In this paper we propose an exact expression for the three-point function (12.11) . The cases of odd 
and even N should be considered separately. If N is an odd number then 

l kl,k2,k ' 3 (ai,a 2 ,a 3 ) if condition (I2.2a[) is true, 
C klMM {a^ a 2 , a 3 ) = { C fcl ' fe2 ' fc3 (a 1; a 2 , a 3 ) if condition fl2T2b|) is true, (2.9) 



otherwise. 



If N is even then both conditions ( I2.2ap and (I2.2bj) can be satisfied at the same time (of course with 



different k and k'). This is possible only if k\ + k 2 + k 3 = 2k > N. In this case our proposal reads as 

C fcl ' fc2 ' fc3 (« 1; « 2 ,«3) = 

' C kl,k2 ' k3 (ai, a 2 , a 3 ) if condition (I2.2ap is true and condition (12.2bj) is wrong, 



C kl ' k ' 2 ' k3 (ai,a 2 ,a 3 ) if condition (12.2bj) is true and condition (I2.2ap is wrong , 
C fcl ' fc2 ' fc3 (ai,a 2 ,a 3 ) + t klMM {a Xl a 2 , a 3 ) if both conditions (l2"2aD and ( l2~2bj) are true, 



v otherwise. 

Functions C kl ' k2 ' k3 (ai, a 2 , a 3 ) and C kl ' k2 ' k3 (ai,a 2 ,a 3 ) in the equations above are defined by (I2.3al) and 
( I2.3bl) respectively. We emphasize that in the definition of C kl ' k2,k3 (ai, a 2 , a 3 ) parameter k is defined 
by 

k 1 + k 2 + k 3 = 2k, 
while in the definition of C kl,k2,k:i (ai, a 2 , a 3 ) by 

h + k 2 + k 3 = N + 2k. 

In particular, in the third line in the r.h.s in (I2.10p C kl,k2 ' k3 (ai, a 2 , a 3 ) is supposed to depend on pa- 
rameter k = {ki + k 2 + k 3 )/2 while C kl ' k2 ' k ' 3 (ai, a 2 , a 3 ) depends on parameter k' = k — y. We note that 
for N = 1 and N = 2 our result reproduces known exact expressions for the three-point correlation 
functions in bosonic Liouville field theory [2H3BII3] and in supersymmetric Liouville field theory (HUTS]. 
One important comment is needed at this point. Throughout this paper we assume that holomorphic 
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and antiholomorphic parafermionic currents commute with each other which is actually not the case. 
The locality of the theory requires 

{z') = uj^{z')^{z), (2.11) 

with u = e 2t7T ^ N . The commutation rule ( 12. lip gives some phases which can be easily restored in our 
final answer ( 12. 9 j) and (12. 10[) . We prefer to hide these phases and put u> = 1 for simplicity. 
Proposed exact expressions ( 12. 9 p and (I2.10p posses a symmetry 

1 

V 

where cosmological constants \i and jx are related via 



b^ T , kj-tN-kj, fi-tp,, (2.12) 



, bQ\ ,2bQ 

% r/i7 ' J N 



7 r/i7 ' ~N~ ' 



Transformation (12.12p reflects the duality of parafermionic Liouville theory (II. ip . Namely, two theories 
one defined by Lagrangian ( II. ip and another defined by 

C = £ PF + *L (d aV f - l~lN ^+^+ e 2b " v , (2.13) 

47T 

are equivalent. Here \l/ + = ^n-i is the conjugated parafermionic current (see appendix [A}. In particu- 
lar, it follows from this equivalence that parafermionic LFT is unitary at b = 1. 

The quantum field theory defined by the Lagrangian density (11.11) is interacting field theory for 
Li 0. However, as was noticed by Goulian and Li [44] in some cases it can be studied by free-field 
methods. By definition, the three-point function ( 12. ip is given by the path integral 

C k ^(a u a 2 ,a 3 ) = [ V^\0)V^\l)V^\oo) e - s ^d^}[d^ (2.14) 



One can perform integration over zero mode ipo of the field <p = (po + ip where field (p has no zero Fourier 
component 

C k ^ k3 ( ai ,a 2 ,a 3 ) = J V^\U)V^\l)V^\oo) (-liN J e 2b ^d 2 z^j e - s ^[dip}[d^], 

(2.15) 

where 

Q — a 

s = — ; — , a = a.\ + a 2 + «3 (2.15a) 
b 

and So[tp, \&] is the action without exponential term i.e. at ll = dE In (12.150 we use the notation 
Va = <Jke 2a ^. According to (I2.15P the three-point correlation function C kl ' k2,k3 (ax, a 2 , 0:3) has a pole 
at s = n (i.e. at a — Q = —nb) with the residue 



-nb TV. 



Res 

ot—Q=—nb 

x f (Vj^O)^!)^ . ..V b {0 \z n ,z- n Mz n Mz- n ))d 2 z 1 . . .d 2 z n (2.16) 



4 Equation (12.14)) needs some comments. First of all the integration over [dtif] is rather schematic because there is no 
known Lagrangian description of parafermionic CFT for N > 2. Second, using projective invariance one can always fix 
the position of any three conformal fields in arbitrary positions [6] (we fix them as usual at 0, 1 and 00.) 

5 The term Q in equation (j2.15aj) comes as usual [4] from the curvature term which is missing in ()1.1[) for shortness. 
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and (• • • ) denotes the path integral with the action without exponential term. The integration J d 2 Zj 
in ( I2.16P goes over two-dimensional plane. For future convenience we represent parameter n in (I2.16P 
as n = mN — k, where m and k are arbitrary integers restricted by m > and < k < N, i.e. we 
define the parameter k as a deficit of the parameter n to be divided by N. We note that correlation 
function in the r.h.s. in fl 2 . 1 6 [) essentially factorizes into bosonic and parafermionic parts. The bosonic 
correlation function in (" 12 . 1 6[) can be evaluated using Wick rules with the pairing 



y?(^Mo,o) 

so that ( 12 . 161) can be rewritten as 



N 



log \z\ + 



Res 



(/iiV) 



mN-k 



(mN-k)b 
mN—k 



n 

i=l 



Zi 



\2A\ 



where 



A 



mN—k 

ir n 

i<j=l 

(2ba 1 + fci 
N 



(mN - k)\ 



X 



ZmN-k) d Z\ . . . (2.17) 



B 



(2ba 2 + fc 2 ) 
N 



and q = — — . 



(2.18) 



(JV|fcl,fc 2 ,fc3). 



z%, . . . , z m N^k) (I2.17P which is actually a polynomial (see below) is defined through 



Function V, 

the correlation function in parafermionic CFT as 



vW k ^\zi,...,z mN „ k ) 



mN—k 



mN-k 



n 



Zi 



Zi 



Zj\ N x 



1|^ [J |; 

i=l i<j=l 

x (a kl (0)a k2 (1)^3(00)^(21)^(^1) . . . *(z mJ v-jfc)t(z mJV -jfc)). (2.19) 



One can easily show (see appendix [A]) that this correlation function is non-zero only if the parameters 
k and kj are related via one of the conditions (12. 2p . Below we will consider in details the case when the 
parameters k and kj are related as in (I2.2al) with the case of the relation ( I2.2bl) being done in complete 
analogy (we present some relevant formulae for this case in the end of appendix EJ). We also denote 
throughout this paper 

n = mN-k, (2.20) 

so in any expressions below parameter n is always given by (I2.20p (except section |4] and appendix [D] 
which have their own notations). Under above assumptions we left with the problem of computation of 
the Coulomb integral 



»(iV[fcl,fe2,&3) 



(A,B\g) 



nl 



x 



n 

8=1 



\2A\ 



\2B 



n 

i<j=i 



,.|4p 



<p(N\kiMM) , 



Zl, 



d 2 z l ...d 2 z n , (2.21) 



which can be viewed as a generalization of the similar Coulomb integrals studied in [1T|B5] in the case 
of bosonic theory (N = 1 in our case) and in [46J in the case of supersymmetric theory (N = 2). It is 
clear that in some domain of the parameters A, B and g this integral is convergent. If one succeeded 
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and computed the integral (12.211) (which is usually given in terms of product of T-functions) then one 
can formally continue it to arbitrary values of the parameter s = (Q — a)/b (see below). This method 
was used in number of papers [4^[TU | [TT j [T4 ^ l4T t l4"g] and in all cases it gives correct expression for the three- 
point correlation function (it is usually confirmed by other methods like semiclassical approximation, 
bootstrap approach etc). 



3 Parafermionic polynomials and Selberg integrals 

The main subject of this section is to compute the integral (I2.2ip . As a warm-up exercise we consider 
contour version of the integral (I2.2ip 

1 i 

» » n n 

3^ k ^\A,B\g) = J ...J Wtfil-U) 8 \{ \t i -t j \^P^ k ^\t ll ...,t n )dt 1 ...dt ni (3.1) 



o o 



i=l i<j=l 

where integers ki, k 2 , k 3 and k are restricted by the condition 

k 1 + k 2 + k 3 = 2k, < kj < k < N, (3.2) 

and n = mN — k. Polynomial P^' fcl ' fc2 ' fc3 )(£ 1? . . . ; t n ) is symmetric polynomial defined through the more 
general polynomial P^' kl ' k2 ' k3 \xi, x 2 , x 3 \ti, . . . ,t n ) as 

Pf 1 ' fe ' fe )(t 1) ..,g = (-l) Mfcl lim a;3 (m - 1)fc3 P^l fcl ' fc2 ^(0,l,a;3|ti,...,t n ). (3.3) 



We stress that the polynomial P^' fcl,fc2 ' fc3 )(i 1; . . . ? £ n ) differs from the polynomial Vm^ kl ' k2 ' k3 ^ (t 1? . . . , £ n ) 
defined by (I2.19P only by normalization 



v^ kM3 \t 1 ,...,t n ) = (-iY m - 1 ^ 



k\(N-k)\ 

— p{h,k 2 ,k 3 ) 



Pff^^l, • • • ,*») (3.4) 



where p{k\, k 2 , k 3 ) can be found in (12.41) . Our choice of normalization will become clear below (see 
appendix IA1 and especially eq (1A.22P ). 

Polynomial P^' kl ' k2 ' k3 \xi, x 2l x 3 \ti, . . . , t n ) is symmetric polynomial in variables ti,...,t n (and not 
necessary symmetric in Xj) satisfying the set of properties which follow from the definition (I2.19P (see 
appendix lAll^l 

P(\xi\\t 3 ) = X^-^Pixiltj), (3.5a) 

3 n 

p (*7% l ) = U x J im ~ 1)kj n^ 2(m ~ 1)p W*;)> ( 3 - 5b ) 

j=l q=l 

P(xi + \\tj + A) = P(xi\tj), (3.5c) 
P(xx,x 2 ,x 3 \x,. . .x,t u . . -tn-JV-i) = 0, (3.5d) 

N+l 

P(x 1 ,x 2 ,x 3 \x j ,...x j ,t 1 , ...tn+kj-N-t) = for j = 1,2,3. (3.5e) 

N+l-kj 



6 In order to save space we will denote Pm fcl,fc2,fe3 ^ X 2, xa\ti, . . . , t n ) simply as P(xi, X2, 2:3^1, . . . , t n ) or even 
P(xi\tj) when needed. 
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Proposition 3.1 Polynomial P^' kl ' k2 ' k3 \x\, x 2 , x 3 \ti, . . . , t n ) is uniquely determined by the set of prop- 
erties (13. 5p up to a total factor. 

The proof of proposition 13.11 is given in appendix [B] The polynomial P^' fcl ' fc2,fc3 ^ x 2 , x 3 \ti, . . . , t n ) 
can be explicitly written in terms of symmetrization of certain basic "monomial" p(xi, X2, x 3 \ti, ■ • ■ , t n ) 
as 

p(N\ k uk2 ,k 3 )( XuX27 X3 \ tu ..., tn ) = A" 1 S\p(x u x 2 , x 3 \t u . . . , t n )], (3.6) 
where S means symmetrization over all variables (ti, . . . ,t n ). Normalization constant 



A, 



N n k\(N -k)\ 



{m\) N - k {{m- l)\) k (3.7) 



in (13. 6p were chosen in such a way that the polynomial explicitly given by (13. 6 p together with (13 .4p gives 
§39}. 

The basic "monomial" p(x\, X2, x 3 \ti, . . . , t n ) in (13 .6p can be defined as follows. We divide our points 
ti, . . . , t n into N groups 

(tx, t-2,-, • • • , t n ) )• (Si, . . . , Sjy), 

where N — k groups contain m different elements and the remaining k groups contain m — 1 different 
elements. We associate with each group s q the "monomial" Yli<jes ~^j) 2 anc ^ take the product over 
all groups 

N 

n n (*-**) a - ( 3 - § ) 

<?=i i<jes q 

It is clear that if the basic "monomial" is proportional to the product (13.81) then the clustering property 
f!3.5dp follows automatically because if one sets any N + 1 variables tj to the same value then at least 
two of them will be inside the same group and hence the expression (13. 8p vanish. The rest of the basic 
monomial is constructed as follows. All groups (si, . . . , s^) are divided into four sets with k — k±, k — k 2 
and k — k 3 groups ((m — 1) elements in each group) and N — k groups (m elements in each group) 

Xl = ( s l> • • • ) s k-ki), Xl = {Sk-ki+1, ■ ■ ■ , Sk 3 ), 
X3 = (Sfcs+l, • • • , Sfe), X4 = (Sfc+i, • • • , S N ). 

We define the basic monomial as 



p(x 1 ,x 2 ,x 3 \t 1 , ...,t n ) = 

N 

= n n & - n & - - n & - x ^ - ^ n & - x ^ - ^ ( 3 - 9 ) 

q=i i<jesq jexi .76x2 jex3 

We note that the expression involving variable Xj in (13.91) looks like 

n^-^-)--- n (^-^)' ( 3 - iQ ) 

i. e. it involves the product over exactly kj groups for any j = 1,2,3 and apparently ensures clustering 
property (I3.5e|) . After symmetrization over (£1, . . . ,t n ) we get the polynomial which is non-zero and it 
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it is easy to see that the rest of properties in (I3.5P are also valid. With the polynomial defined by ( 13. 6p 
and ( 13.9P we define the polynomial Pm^ kl ' k2 ' k3 \ti, . . . ,t n ) as in ( 13.31) . From the standard CFT point of 
view the definition ( 13. 3p is equivalent to using the projective invariance of correlation function 

{(Jk x {xi)(JkA X 2W 3 {x^{tx) . . . -> («7 fcl (0)t7 fca (l)<7 fcs (oo)*(ti) . . . *(t n )) s 

allowing to fix x x — 0, x 2 — 1 and x 3 = oo. 
Integrals of the form (13.11) i.e. 

i i 

» » n n 

j ... j ntfii-tjf n it i -t i i 2 «p(t 1J ... J t n )tzt 1 ...dt n , (3.H) 







j=l i<j=l 



where P(t x , . . . ,t n ) is some symmetric polynomial were extensively studied in the literature. When 
P(t x , . . . ,t n ) = 1 the integral ( 13.111) coincides with ordinary Selberg integral [12]. The case when 
P(t x ,...,t n ) is elementary symmetric polynomial was considered by Aomoto [30] • Further gener- 
alization has been computed by Kadell [50]. In this case P(t x , . . . ,t n ) coincides with Jack polyno- 
mial Jac]/ 9 (ti, . . . where A is arbitrary partition of at most n parts (for further development see 
also [51]|52]). In this paper we present another exact results for the integrals of the form (I3.1ip (see 
below). 

We will follow mainly original approach of Selberg [32] also described in [53]. For g being non-negative 
integer number we can expand in (13.11) 

n 

n i** - p^^ih, . . ., tn ) =j2 c ^ # ■ ■ -c, (3-i2) 

i<j = l Uj 

where the sum goes over all eligible sets of integers . . . ,u n . It is clear that the values of the coefficients 
Cvi,...,u„ in (I3.12p are independent on the order of the integers V\, . . . , u n . It happens that in complete 
analogy with Selberg case (i.e. in the case with polynomial being a constant) the integers arranged as 
V\ < Vi • ■ ■ < v n satisfy proper set of inequalities which allow to find A and B dependent part of the 
integral (13. ip analytically. 

Proposition 3.2 For V\ < ■ ■ ■ < v n the non-zero values of C Ul Vn in (13.121) occur when the integers 
VpN+j satisfy the inequalities 

VpN+j > P+ (pN + j - 1)5- if p = 0, ...,m-2, j = 1, . . . , JV - k x , 

u pN+j >p + l + { P N + j-l)g if p = 0,...,m-2, j = N - k x + 1, . . . , N, (3.13a) 

VN( m -i)+j > m- 1 + ((m- l)N + j - l)g if j = I, . . . , N — k, 

VpN+j < m + p - 1 + ((m + p)N - k + j - 2)g if p = 0, . . . , m — 1, j = 1, . . . , N — k, 

VpN+j <m + p-l + ((m + p)N -k+j-2)g if p = 0, . . . ,m — 2, j = N — k + 1, . . . , N — k + k 3 , 

u pN+j <m+p+({m+p)N -k + j -2)g if p = 0, . . . ,m - 2, j = N - k + k 3 + 1, . . . , N. 

(3.13b) 

We give a proof of proposition 13.21 in appendix O Substituting expansion ( 13.121) in the definition ( 13. lj) 
of the integral j£ lkl ' k2 ' ks) (A, B\g) and using the Euler beta integral formula we get 

jjjn (A, B\ g ) = JJ c, „„ f[ ^na + t + B^.f' - < 3 - 14 ) 

Vj j=i ^ ^ ^ V 
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Using the first set of inequalities (13.13aj) in the numerator in (I3.14p the second set (13.13bl) in the 
denominator and the relation T(l + x) — xT(x) we can rewrite 

TT ni + A + v,) _ R yr 2 G^\l + A + P Ng + p) 

IIVQ + A + B + Uj) { ^ }V ' l} QG fJ k3{1 + A + B + ((m+p)N-k-l)g + m + p) 

x G^- k) (m + A+(m-l)Ng) 

G ( N ~ k} {2m + A + B+{{2m-l)N-k-l) g y 



where 

N-k N 

Gi N \x)^Y[T(x + (j-l)g) J] r(l + x + (j-l)g), (3.16) 

j=l j=N-k+l 

and H{ Vj y(A, B) is some polynomial in A, B and g (which depends on the integers {vj} = (ui, . . . , u n )). 
By symmetry reasons it is also useful to represent in (13.141) 

n 1 m— 2 

H r(l + B) = —— H G[ N 2 \l + B + pNg + p) G { N - k) (m + B + (m - l)Ng), (3.17) 

j=l yk 2 {&) p=Q 

with some polynomial Qk 2 {B) (it does not depend on choice of integers {vj}, but its explicit form is not 
important for us now). Using (13. 15ft and (13. 17ft we get 

j(N\mMto) (A B\a) = H{AB) TT G ^ 1 + A + P N 9 + P) Gi £ ) ( 1 + B +P N 9+P) x 

Qk 2 (B) LL G f) k ^i + A + B + ((m + p)N-k-l)g + m + p) 

x Gj N ~ k \m + A + (m- l)Ng) G^jm + B + (m - ljiVg) 
GQ N ~ k \2m + A + B + ((2m — l)iV — k — l)g) 



where 

Now due to the symmetry of the integral (13. ip with respect to A -H- 13 and fci <(->■ &2 we can also represent 
it as 

j( N lkl MM) (A m = EiAB) yv G^\l + A + pNg + p)Gl i :\l + B + pN g+ p) 

QkM) j;} Gi N X(l + A + B + ((m + p)N-k-l)g + m + p) 

x G { N - k \m + A + (m- l)Ng) G^%n + B + (m - ljjVg) 
G { Q~ k) {2m + A + 5 + ((2m - l)iV - k - 1)#) 

with some another polynomials H(A,B) and Qfci(^)- The validity of both equations ( I3.18P and (I3.19P 
requires 

H(A, B) _ H(A, B) 
Q k2 {B) ~ Q kl {A) ' 
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which means in particular that H(A,B) is completely divisible by Qk 2 {B) as well as H(A,B) is com- 
pletely divisible by Q^A) 

H(A,B) _ H(A,B) 

where R(A, B) is some polynomial in A and B of certain degree^. Let us show, that in fact R(A, B) 
does not depend on parameters A and B. First, we consider large A asymptotic of the integral (13. ip . 
It can be done by changing variables = e~~z~ in (13. ip 

jWkiM,k a )^g^ ^ A-n(m+B+(n-l)g)+(m-l)(k-k 2 ) ^ (3.20) 

From other side using 

— — — — ~ a; , as x — >■ oo, argfa;) 7^ 7 r ) 

in ( I3.18P we get 

G^l+A + pNg+p) 



n 



£ GX\ (l + A + B + ((m + p)N-k-l)g + m + p) 

G ( N ~ k) {m + A + {m-l)Ng) 



(N) 

3 

(N-k), 

^-n(m+B+(n-l)g)+(m-l)(k-k 2 ) /g 

" Gf~ k \2m + A + B + ((2m-l)N-k-l)g) 

Because both (I3.20p and ( I3.2ip have the same behavior at A — > oo we conclude that the polynomial 
R(A, B) does not depend on A. Analogously we find that it does not depend on B as well. So R(A, B) 
may depend only on parameter g, i.e. 



i}oG ( ill 3 (l + A + B + ((m + p)N-k-l)g + m + p) 

x gg^V + A + (m - l)ATg) G^- fe) (m + g + (m - 1)^) 
G { N ~ k) (2m + A + B + ((2m - 1)N - k - l)g) 

The factor Cm ,k \g) in ( 13.221) depends on parameters ki, ki and k% only through their sum i.e. it 
depends only on parameter fc. This fact can be easily seen if we take 

1 Sk . L 2 , 

t k = t, A = B = —, 3.23 

2 2L 2 y J 

consider the limit L — > oo in the integral ( I3.ip ^l and from other side use Stirling formula in the r.h.s. in 
(I3T22D . One obtains that C^' k) (g) is given by generalized Mehta integral 

oo oo 

—J J... J l[e-^l[\s i -s j \ 2 ^^( Sl ,...,s n )ds 1 ...ds n , (3.24) 



j=l i<j 



7 It may depend on parameter g in non-polynomial way, which is actually the case (see below). 
8 While performing the limit of the integral (|3.ip we used the following limit 



lim ( 1 — ( — 1 | = e x 2 



L— >oo V \ L 
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here Vfim (si, . . . , s n ) is the polynomial given by symmetrization of the "monomial" p(si, . . . , s n ) 

q3^(si, • • • , s n ) = A" 1 S [p(s u s n )} , (3.25) 
where normalization factor A m is given by ( 13. 7p and 

N—k-l m AT-l m-1 

p(si, . . . , S n ) = J^J (Si+mp — Sj+ mp ) J~J J^J {Si+(m~l)p — Sj+( m -l)p) . 

p=0 i<j=l p=N—ki<j=l 

We see that (I3.24p indeed depends only on parameter k and not on parameters kj. Polynomial (I3.25P 
is symmetric homogeneous polynomial of total degree (m — l)(n — k) (we remind that n = mN — k) 
satisfying clustering property 

^(x 1 _^x / , s u . . . , s^.x) = 0. (3.26) 

N+l 

We note that the polynomial (13.251) can be obtained from the polynomial Pm^ kl ' k2 ' k3 \ti, . . . ,t n ) with 
variables tj taken as in (I3.23P in the limit L — > oo. As was shown in [SUES] the degree (m — l)(n — k) 
is the minimal possible degree for the symmetric homogeneous polynomial in n = mN — k variables 
satisfying (13.261) and this unique polynomial is proportional to Jack polynomial [56] 

¥£\s 1 ,...,s n )~J ac a x "(s 1 ,...,s n ), (3.27) 

where = —N — 1 and A is the special partition 



A 



2(m - 1), ... , 2(m - 1), 2(m - 2), ... , 2(m - 2), . . . , 2, . . . , 2 

S v ' V v ' S v ' 

N-k N N 



Since the factor Cm' k \g) does not depend on parameters kj we found convenient to evaluate it from 
the integral ( 13. ip in the case k\ = and k 2 = k 3 = k. By Carlson's theorem [570 we find that « 
holds not necessary for integer values of the parameter g. In particular it can be continued to the region 
of sufficiently small negative g. In this region function Cm \g) can be found if we change variables 

tj=T^ With J2^ = 1 fOT J = 1 --- N i 

in the integral (13.11) with k% = and ki = k% = k 

N 1 1 N 

J jd^...j d^f[^{l - i 3 r) B J] % - ^S(^ + • • • + - !)• 



j=l i<j 
n—N N n-N 





1 1 



/n IL 1 V 1 V IL 1 V 

... / \[tf{l-U) B \[\U-T^ J] \t i -t^P^ k \r^...T^ N ,t U ...,t n . N )dt 1 ...dt n . 

o i=l i =1 i <i= 1 

(3.28) 



9 The only condition we need is that r.h.s. of (|3.22[) grows as 0(e M ' 9 ') with /i < %. This assumption will be justified by 
explicit form of Cm (g)- 
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The integral over dr in (I3.28j) is divergent logarithmically at A = —1 — (N — l)g. Applying the general 
formula 

lim (1 + a) / x a f{x)dx = /(0), 

Jo 

for any function f(x) continuous at x — > + we can consider the limit 



lim (l + A + (N-l)g)<MB. 



Using the clustering property 



P^ix^tu . . .,t n _ N ) = fwU (tj - xypw^h, . . .,t n _ N ), 
n 2 i =1 

which follows from the explicit form of the polynomial Pjf' ' (ii, . . . ,t n ) one immediately obtains 
lim (l+A + {N-l)g)j£^\A,B\g) = N-% r ^^ 

A-t-l-(N-l)g {n — i\)l 

(3.29) 

where 

i i N 

X N (g) = i I . . . I f[ tf-W-Di JJ (6 _ {.psfa + . . . + Cjv _ i) ^ . . . rf6v . 

o o fc=1 *<J 

In order to compute Xjsf(g) we set in (13.291) m = 1 and = 0. In this case the polynomial entering the 
integral °' ' ^ (A, -B|g) is actually some constant so that J ( f f ^ ' '°\A, B\g) is ordinary Selberg integral 
which is already known. Considering the limit (13.291) in this case we get 

Xn{9) - n^i) l\ m+g) • (3 - 31) 



Now substituting (I3T22D in f[3~29j) we find that 

,(AT) 



r (iV,t]M _ Ar-f " ! G iv-fc(™ -l + ((m-l)N-k + l)g) r ( N ,k) ( , 



(3.32) 



Repeating f l3.32|) (in — 2) times we get 

im-p* n\ j^ 1 G$l k (p +(pN-k + l)g) 
p=i 

Finally using that under chosen normalization 



, fJV!(iV-*)! 



and hence the integral (I3.24p is just the ordinary Mehta integral in this case we find that 
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Combining altogether one arrives at 

Cg' fc) (g) = n! [^- w fc , (iV l fc)! J (faT^yj G ° + ^ l[Gfl k (p+(pN-k + l)g), (3.35) 

where n = mN — k. Finally we conclude that the integral ( 13.1 j) is given by fl3T^|) with C^f' fc) (</) given 
by (UBS)). 

Now we are ready to compute the integral Qm^ kl ' k2 ' ka \A, B\g) defined by (12.21 j) . One can use the 
following observation. For any integral over the plane (more precisely over n two-dimensional planes) 

[■■ III \M 2 % - 1| 2S II I** - *il 49 i P (*i> • ■ • > f n)\ 2 d 2 ti ■ ■ ■ d2t n, (3-36) 

J J i=l i<j 

where P(ii, . . . ,t n ) is arbitrary polynomial the following relation is hold [IT j lo'S j lo'U] 

I n = X n {B\g) lj n , (3.37) 



I, 



where 



pi i-i n 

I n = l\tf{l-t i ) B T[\t i -t j \ 29 P{t 1 ,...,t n )dt l ...dt n , (3.37a) 

Jo J ° i=l i<j 



and 



/OO /'OO n 

■ ■■ / Y[tf{t i -i) B Y[\t i -t j \ 29 P{t 1 ,...,t n )dt 1 ...dt n 
^ 1 1=1 i<j 



(3.37b) 



In fact (I3.37P was proven for pure Selberg case (i.e. with P(ti, . . . ,t n ) = 1) by standard manipulations 
with contour integrals [60J. Since the polynomial is a single- valued function it is not an obstruction to 
deform contours in a same way as was done for P(t\, ... ,t n ) = 1. Moreover, since the factor X n (B\g) 
originates from the non-analyticity of the integrand in (I3.36j) it is going to be the same for any polynomial 
P(t 1 , . . . ,t n ). Explicitly X n (B\g) reads as 

X n (B\9) = ff sinn{B + j9) ^l+M . (3.38) 

n\ - LJ - sm7rq 

In order to compute the integral (I3.36|) one has to compute the contour integrals f)3.37al) and (I3.37bj) . 
The late can be reduced by change of variables tj — > 1/tj to the integral of the form (I3.37al) with some 
another polynomial P(ti, . . . , t n ) defined as 

n 

P(t l ,...,t n ) = l[t]P(l/t 1 ,...,l/t n ), 

i=i 

where 77 is some integer number. In the case considered in this paper polynomials P(ti, . . . ,t n ) and 
P(ti, . . . , t n ) belong to the same class. Namely, as follows from f l3.5bl) 

n 

Ilf^ P ^ lkM \l/h, l/t n ) = P^ k ^\t u . . .,t n ). (3.39) 

3=1 

16 



Using (I3.37p . (13.38j) and ( I3.39j) one can compute the integral (12.211) 

m— 2 

X S 7^ (p** + ^P^) ^ + 7 g0 (2^ + ^2) 7 W ( 2^ + ^) X 

1 

7? (f) 7r _fc) + s - 1)^ 2 ) 7f + (m - 1)& 2 ) 7r- fe) + (™ - 1)^) 

(3.40) 

where integers ki, k 2 , &3 and A; satisfy the condition (I2.2a|) . In ( 13.40}) the factor p(ki, k 2 , k 3 ) is defined 
by (12.41) . We note that the parameters A, B and g are expressed in terms of the parameters atj and b 
as in (12. 18}) while the parameters act, a 2 and a 3 are supposed to be related by the screening condition 
(we recall that n = mN — k) 

q;i + a 2 + «3 + nb = Q. (3-41) 
In (13.401) we used also the following notation 

7i^n\(*+^^) n (3 - 42) 

3=1 ^ ' j=N-k+l ^ ' 

and 7^ k \x) = Yif=i "f(x + (j — l)bQ/N). Now it is easy to show that the function C kuk2,k3 (ai, a 2 , a 3 ) 



defined by eq (I2.3aj) satisfies 

Res 



C kM3 { ai ,a 2 ,a 3 ) = g^ hM *\A,B\g), (3.43) 

a— Q=— nb 



where Q^ lklM ' k3 \A, B\g) is given by (CT) . 

In the dual case when parameters k\ , k 2 and are related as in (12.2bj) corresponding integral can be 
calculated in a same way (we present some relevant formulae in the end of appendix |A]). The resulting 
integral satisfies the residue condition similar to (13.431) with function C* 1,fc2 ' 3 (oti, a 2 , a 3 ) given by ( I2.3bl) . 
We note that if N is an odd number then only one condition (12.2a)) or (12.2b)) can be satisfied for given 
k\, k 2 and k$ and hence the three-point correlation function (12. ip (if it is non-zero) is given ether by 
(I2.3al) or by (12 .3bj) as in (12. 9p . Contrary, for even N both conditions (I2.2al) and (I2.2bj) can be valid at 
the same time. In this case it is natural to expect that the three-point function will be given by sum 
of (12.3a)) and (12.3b)) . The phenomenon of this type appears already for N = 2 i.e. for supersymmetric 
LFT for the three-point function of two Ramond and one Neveu-Schwarz fields [15]. Our proposal for 
even N is given by (12 . 10)) 

We note that the residue condition (13.43)) and similar condition for the function (12.3bl) are not enough 
to determine the three-point function C kl > 2 > 3 (a\, a 2 , a$) completely. Moreover, as one can easily see the 
proposed exact expression (12.9l) - (12.10p has much more poles than predicted by ( 12. 15ft . Our conjecture 
is inspired by intuition adapted from the bosonic [lj and supersymmetric [14J Liouville field theories 
and has to be justified by independent methods. In the next section we perform non-trivial check of 
our conjecture based on dual representation of the parafermionic LFT by the three-exponential model. 
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4 Dual representation of parafermionic LFT 



In this section we consider dual representation of parafermionic LFT based on three-field model sug- 
gested in [39]. Some notations which we are going to use in this section will overlap with notations used 
before. 

We consider the theory with the action, which contains three exponents 



S + _ = M exp{ai(pi - a 2 y?2 - ia^z}, S_ + = Mexp{-aiy?i + a 2 f2 - ia^s}, 
#++ = Mexp{a 1 (p 1 + a 2 f 2 + iaip 3 } 



(4.1) 



where parameters ai,a 2 and a satisfy the condition: 

a 2 - a\ - a\ = 1/2. (4.2) 



Namely 



3 



The affine version of the theory (14. 3 p where the term is added to the action was considered in 
different regimes^ in [32JEI1E2]- The theory ( 14. 3 p is non-unitarian CFT with the stress energy tensor 



and central charge c = 3 — + ^ + -^z. It is convenient to write our parameters a, otj in the form 



W = N + 2, 4a? 4a 2 2 = T fp- (4-5) 



Then the central charge will have a form 



3N + (4.6) 



N + 2 N 

which is exactly the same as (11.21) . The primary fields will be exponential operators 

V(a) = exp (oi^i + a 2 ip 2 + ia 3 tp 3 ) (4.7) 

with conformal dimensions 

A(a) =*iU-~ «x) +^(^r- «*) - «a - a 3 



10 There are three regimes in the affine theory when all parameters a and otj are imaginary, when a is real and both otj 
are imaginary and finally the regime considered in this paper when all parameters are real. In the first regime when all 
parameters a and ctj are imaginary the theory is unitary and its spectrum consists of quadruplet of fundamental particles 
and their bound states |39j . The theory is integrable with two-particle S- matrix being proportional to the tensor product 
of two-soliton ^-matrices of the sine-Gordon model. The integral representations for the form-factors of local fields in 
this theory were studied in [61,. The particle content and scattering theory in the remaining non- unitary regimes of affine 
theory are rather complicated and deserves further studies. 
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The correlation functions of exponential fields can be calculated using Goulian-Lie method described 
above 

nrf-\t \ \ r ( ~* \ ( w r(n 2 - m - n 3 ) T(ni - n 2 - n 3 ) T(ni + n 2 - n 3 ) 
(V(a 1 )(x 1 )...V(a r )(x r )) = -± x 

OL\ a 2 a 

x <F(aW)(a: 1 )..^(aM)(x r )(-S;_)-(-S: + )™ 2 (-S+ + )-> FF (4.8) 

where S+_ = M J cPxS+_(x)... and index FF denotes that expectation values are taken over free field 
vacuum. The numbers are restricted by the conditions 

r 1 r 1 

S~] a[ 3) + (m -n 2 + n 3 )ai = - — , V] a 2 j) + (n 2 - ni + rc 3 )ati = - — , 
^-f 2ai ^— ' 2a 2 



a 3 + (ra 3 — ni — n 2 )a = — . 



1 

2 a 



Free-field correlation function in the r.h.s. in fl4.8[) can be expressed in terms of Coulomb integral only if 
all the parameters are non-negative integers otherwise eq (14.81) is rather formal. The calculation of the 
correlation functions in our theory is very similar to the same calculation for sine-Liouville model [64] 
(see also [55]) and based on usage of certain integral identity (see for example eq (1.13) in [47]). We 
prefer to hide details here and present only final answe Let us suppose the condition 

E«?-^ = r'a (4.9) 
i=i 

with some integer number r' and define p as 

p = r-r'-2, (4.10) 

then correlation function of the fields V(dj)(xj) can be expressed in terms of correlation functions of two 
Liouville field theories with coupling constants 2ct\ and 2a 2 with cosmological constants — /i 2 7(— 4a 2 ) and 
— /i 2 7(— 4a 2 .) (we label corresponding correlation functions with lower subscripts 1 and 2 respectively). 
Namely, 

(V(a 1 )(x 1 )...V(a r )(x r )) = ^ f[N{aj) J d% . . . d% Q(t u . . . , t p \ Xl , . . . , x r ) 

(v. ai (ti) . . . v„ ai (t p )v a ( 1)+ai (^i) • • • v^^XrMV-^h) . . . v_ a2 (t p )v a w +a2 ( Xl ) . . . v a ir )+a2 (x r )) 2 

(4.11) 

where V ai are the standard exponential operators in Liouville field theory with conformal dimensions 
a>i(Q — CL-i) and factors Af(aj) are 

Af(dj) = 7 (1 + 2aa ( 3 j) - 2 ai a[ j) - 2a 2 a ( 2 j) ). (4.12) 



interested reader will be able to obtain (|4.11l) with the help of integral identity (1.13) from [47] in a simple way. 
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The factor Q(t x , . . . , t p \x%, . . . , x r ) in equation (14. lip is given by 

Q(t 1 ,...,t p \x 1 ,...,x r ) = — p r fc<m . (4.13) 

nni^-^i M4m)+Q) 

i m 

We see that our model, similar to sine-Liouville model admits the separation of variables and possesses 
the dual sigma model representation (see appendix [D]). The general three-point correlation function in 
our CFT is rather complicated object which we suppose to study in separate publication. However it 
is rather easy to calculate the two point functions which coincide with special reflection coefficients of 
the fields V(a). We normalize our fields V(a) by the condition 

V(ai, a 2 , a 3 )V - a u - a 2 , ^ - a 3 j ^ = 1 (4.14) 

Then the reflection amplitude defined as V(ai,a 2 ,a 3 ) = R\(a)V — ai , «2 5 coincides with the 
two point function 

Ri(ai, a 2 , a 3 ) = (v(ai, a 2 , a 3 )V (on, ^ a 2 , - a 3 ^j \ . (4.15) 

In the similar way can be defined amplitude R 2 (a) and amplitude associated with double reflection 

V(ai, a 2 , a 3 ) = i?i j2 (ai, a 2) a 3 )V ^— ^ a x , — ^ a 2 , a 3 

which coincides with the two-point function 

i2i )2 (ai, a 2 , a 3 ) = (v(ai, a 2 , a 3 )V (ai, a 2 , ^- - a 3 jj (4.16) 
The amplitude associated with full reflection 

/ 1 1 i 

V(ai, a 2 , a 3 ) = R(a)V - a u - a 2 , - a 3 

\2a% 2a 2 2a 

is equal to the two-point function 

R(ai, a 2 , 03) = (V(a u a 2 , a 3 )V (a x , a 2 , a 3 )) (4.17) 
It is easy to derive from eqs. (j4.11|4.12p that functions Ri, R 2 and Ri t2 have the form: 

R\(a) = — 2ai \ , — r, 4.18 

7(2«iai + 2a 2 a 2 — 2aa 3 )^{2aiai — 2a 2 a 2 + 2aa 3 ) 

R 2 (a) = - gy^gj (419) 

7(2«iai + 2a 2 a 2 — 2aa 3 )'j(— 2a x a x + 2a 2 a 2 + 2aa 3 ) 
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and 

Rl 2 (a) = JCK+ °.)*U* + °*) (4.2,,) 

7(2aiai + 2a2a.2 — 2aa3j7(— 1 + 2ajai + 2a2fl2 + 2aa3j 

Here the functions i?^, are the reflection amplitudes (two-point correlation functions in the Liouville 
CFT with coupling constants 2a j and cosmological constants /i, = — 7rM 2 7(— 4a?)) 

/vrM\ { ^" 4ai)/2Ql r(4a 4 (ai - tMIYI + i(a< - ^) 

*S*(* + <*) = S FTT^ ¥^ 4a lV (4.21) 



It is easy to see that reflection amplitudes satisfy the "associativity condition" 

#1,2(^1, Oi, a 3 ) = -Ri(ai, a 2 , a 3 )i? 2 (7^ a x , a 2 , a 3 ) 

zai 

The reflection amplitude R3 defined as V(ai, a 2 , 03) = i?3(a)V (ai, 02, j- — 03) coincides with correlation 
function 

i2 3 (a) = /v(a 1( a 2 , a 3 )^ ( — «i, 7^ «2, 03 

It can be calculated with the result 

nM\ ^ )/2a r(-4a(a 3 - £))r(l + l -{a, - £)) 



^ ^ U« 2 ) r(4a(a 3 - i))r(l - i(a 3 - £)) X 

x 7(1 — a 3 a 3 + a 2 a 2 — aa 3 )7(l — a 3 a 3 — a 2 «2 + aa 3 ). 
The amplitude R (a) = _R 3 (a) i?i,2 (a)- If we introduce the momentums p 1( p 2 and q by the relation 

2aj-aj = - + i2ajPj, 2aa 3 = — + 2ag, 

then denominator in eq.f l4.20p will have a form 

7(- + 2za!P! + 2za 2 p 2 - 2ag)7(^ + 2ia 1 p 1 + 2ia 2 p 2 + 2aq), 

and it easy to see that amplitude -Ri^ will be unitarian with respect to pj —> —pj. 

At this point we will not discuss here any more the general CFT (14. 3 ft and start to consider the 
special subspace of our model which is related with parafermionic Liouville model considered above. 
We have seen already that the central charge of our model (14. 6 p coincides with the central charge of 
parafermionic Liouville model (11. 2p . We show that certain three-point correlation functions (12. ip in 
parafermionic LFT can be expressed in terms of special correlation functions of our model. It is well 
known that the parafermionic algebra admits the free field representation in terms of two fields <pi and 
02 with the stress energy tensor [66] 

Tp F (^) = -UdAif ~ ](d z <P 2 ? + ^d 2 z <p 2 . (4.22) 
4 4 4a 

Parafermionic currents ^(z) and ty + (z) can be represented in terms of the chiral parts of these fields 
as [66] 

,'N + 2\ 1/2 ( n ( N N 1/2 



iy j y' d ^+{wT2j ^i)aq>(-^iW/V^), 

1/2 / 1/2 \ ( 4 ' 23 ) 
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The order fields a k can can be represented in terms of these fields as 

-=M-SK-^)> (4 - 24) 

where is the normalization factor which will be specified below. To establish relation between two 
models it is important to express fields (pi, ip from (14. 3[) in terms of the fields ipt from ( 14 . 2 2 [) and vice 
versa. For that it is convenient to compare the stress energy tensor (I4.4j) with stress energy tensor of 
parafermionic Liouville model 

T PFL = T PF - N(dip) 2 + Qd 2 tp (4.25) 
It is easy to see that field (p 2 can be identified with field <^3- For other fields we derive 



2«i0i 2a 2 0i 
tpi = 4a 2 (f (f2 = Aomp + 



N " " y/N ' 



(4.26) 



N N VN VN 

It is interesting to express the screening charges in parafermionic Liouville model which are 

Ji = ^(z)e 2b ^ z) , J 2 = ^ + (z)e 2b ^ z) , 

in terms of fields We derive using the relation 6 = ^ that 



= ( T^) ( iad z^>3 + aidzVi ~ a 2 d z (p 2 ) exp ( rl " > 



AN J y ' ro ' 1 \ 2a 2 

1 \ ( { z) 

J2= ( J (i°tdzV3 - tti^^i + o£ 2 d z (p 2 ) exp f ^ 



(4.27) 



It is easy to check that fields Ji and J 2 commute with all screening charges S+_... given by eq ()4.ip 
and operators j (Px^JiJi can be used as additional screening charges for calculation of the correlation 
functions. We note that there is one additional screening field J3 

J3 = (jj^) (-iad z (f 3 + aid z cpi - a 2 d z (p 2 ) exp ^ • ( 4 - 28 ) 

Now we can consider the simplest three point correlation function in parafermionic Liouville model. 
Namely, C kl,k2,ks (ai, a 2 , 03) where ki + k 2 + k 3 = N and take for simplicity iV odd. In this case r' = — 1 
and hence p = so our correlation function does not contain any integrals. It follows from eqs (I4.26j) 
and that the field V a (fc) has a form 

y , k) = Nfcexp ^_k(a m - a2 <p 2 )^ ^ ^(2^ + 2^ 

The normalization factors Af(a) in eq (14. lip in this case are 



7((8aia 2 aj + 4fcja|)/iV) 
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and the there-point correlation function C kl > k2,k3 (ai, a 2 , a 3 ) with k\ + k 2 + k 3 = N can be written in the 
form 

3 

C kM (a u a 2 ,a 3 ) = J^N^^K) x C 2ai {A u A 2 , A z )C 2a2 {B x , B 2 , B 3 ) (4.29) 

8=1 

where C 2oii denotes the three-point correlation function in Liouville CFT with the coupling constant 
2«j and 

_ 2a 2 a { + a x (N - h) _ 2a x ai + a 2 (N + ki) 



N N 
One can use the following easily established relations 

kt 1 w { 4a 2 a + 2a x k\ (Aa x a + 2a 2 (N + k) 
■N k {a)T 2ai 77 T 2a2 



N / V N 



'Aa 2 a + 2a 1 (N- k)\ (Aa l a + 2a 2 k\ 
KT 2ai 77 i 2Q2 (4.30) 



AT / — 2 \ N 



2a 2 a + 2a 1 k _ a 2 +l\ ( 2a x a + 2a 2 (iV-fe) _ a\ + 1 



and 

T^(a,fc) = (2a 2 ) V - - ) (2a x 

x T I , 2ai 1 T I , 2a 2 ) , (4.31) 

2 

where k = (2a 2 ) 77 , a x and a 2 are taken as in (14. 5p . (x, b) is given by (12.61) and T(x, b) is 

ordinary Upsilon-function. Using (I4.30p and (14.31 j) we find that if we take into account that cosmological 
constants are equal to — , y(—Aaf)M 2 as well as the relation between the parameters fi in parafermionic 
LFT and M in three-exponential model derived in 



l 2 . 

where ^ = and also the explicit form of normalization factors N*. 

'N-k + 1 



iV + 2 



computed in [52] we obtain that (I4.29P coincides exactly with the expression (12.3bp for the three-point 
correlation function proposed above in the case k\ + k 2 + k 3 = N. 

For the calculation in the more general case when parameters kj are related by the condition k\ + 
k 2 + k% = N + 2k we should introduce k additional screening charges generated by the field J3 (I4.28p . 
The calculations become more involved but they give exactly the same answer as in f!2.3bp . 

Now we can calculate the reflection amplitude for parafermionic Liouville field theory. The simple 
calculation gives a result that function defined by the relation 

= Rk(a)V<?_- a k) (4.33) 
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is equal 



These functions satisfy the condition Rk(a)RN-k(Q — a) = 1. For arbitrary b functions Ra(a) = -Rat (a) 
and Rn_(cl) (for even N) satisfy unitarity conditions Ro(a)Ro(Q — a) = i?jv (a)-Rjv (Q — a) = 1. For 6=1 
all functions -Rfc(a) satisfy this condition and full theory is unitary Using this reflection amplitude we 
can derive all three-point functions from the three-point function with k\ + ki + k$ = N + 2k. In this 
way we have an alternative derivation of the result obtained in section [3] 

5 Concluding remarks and open problems 

In this section we discuss briefly some subjects and questions which were not considered in the main 
body of the paper. 

1. The important question in the analysis of any CFT is the description of the chiral algebra of 
theory. Besides the non-local symmetry generated by the field ( ID. Ill) our theory possesses the 
infinite dimensional local symmetry generated by holomorphic W— current W& and other currents 
with even spin which appear in the operator product expansion of this current with itself and 
with these new holomorphic currents- 12 !. These currents form the closed W algebra. For three 
exponential model considered in Section 4 the explicit form of these currents in terms of three 
fields can be found in [521E2JE3]- The fields 

Vi Kq) = *f M) e 2ov (5-1) 

where fields ^ are defined in appendix A are the primary fields of this W— algebra. The full 
space of fields in parafermionic LFT, local with respect the fields V£ can be derived by application 
of this W algebra to the primary fields 

N-l N-l 
fc=0 q=0 

2. In this paper we calculated the three-point correlation functions for the fields V*. Three point 
correlation functions of more general local fields Va k ' g ' > = ^\ qq )& 2aip 

{V^\Q)V^\l)V^\oo)) 

can be calculated by generalization of our method for this case. The analog of the polynomials 
p(A r |fei,fc 2 fc3) f or thjg case can derived by the proper fusions of the points in this polynomials 
dictated by the parafermionic algebra. Namely any field q ^ can be obtained in the operator 
product expansion of "order" field with parafermionic currents \& and ^ '. As an example we 

12 The field W4 and other W— currents appear in the OPE of two fields (|D.11[) in the sector corresponding to identity 
operator. 
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give here the three-point correlation functions of two operators V$ with operator Va k ' q ^ which 
determine the operator product expansion of two fields Va^ (with k\ > k 2 and k\ + k 2 < N) 

k\+k 2 N—l 
a j=ki—k2 g=0 

The structure constant C^'g 1 ' 02 *^ in this expansion can be expressed through the three-point 
correlation function (V^ 1 (ti)V^(l)V^ k+2q) (oo)) . For h + k 2 + k = 21 < N it has a form 



(v£(0)V k >(l)V^ k+2q \oo)) = A(a,k,q,b) 



I bQ\ ,_26Q 



p(k 1 ,k 2 ,k)x 



Tr ) WTir ) (2a 1 )Tif(2a 2 )Ti^(2a 3 ) 

X 



T^(«! + « 2 + a 3 - g)T[ A 2 i+g (a 2 + a 3 - a^T^^ai + a 3 - a 2 )T, ( 3_,(ai + a 2 - a 3 ) ' 

(5.4) 

here function A(a, k, q, b) for < q < N — k has a form 

(jV-fc- g )!jV!* (-JV^ <2,<iV-A; 

g!(AT-A;)!(Ar-g)!2 l = l (2a + {N - k - 2q)b + jQ) 2 ' q ~ 

and for JV - k < 2q < 2N - 2k 

(N-k-q)\N\ 2 N ^r q (-N 2 Y 

~ n\(N - 1AUN - n\W 11 



q\(N - k)\(N - q)\ 2 A ^ (2a + (N — k — 2q)b + (N — k — j — 1)Q) 



2 ' 



For q > N — k the function A(a, k, q, b) can be calculated using the "duality" relation 

A(a, k, q, b) = A(a, N - k,N - q,\). 

b 

3. The important role in studying of parafermionic LFT play the degenerate fields V^ nb , 2 _ m , 2b wnere 
non- negative integers n,m satisfy the relation m — n = k (mod TV). The OPE of fields V* (and 
Va k ' 9 ^) with these fields give only finite number of field V^ k ' 9 \ The four point correlation functions 
of the degenerate field with three arbitrary fields can be represented by the finite dimensional 
integrals and satisfy the differential equations. As an example we give here the simplest four point 
functions in parafermionic LFT. Let i + j + k = I where 21 < N. Then 

<<(oo)v£(i)v£^ 

where 5, = { -^(2a lb + $±$), 5 2 = ^(2a 2 b + §±f ) and 

N-l N-l 



'N-l 



p=l m<n 
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with 



b(a-2a 1 -Q -k/b- {N-l)b/2)_ _ b{a - 2a 2 - Q - j/b - (N - l)b/2) 

N 3 N 

b(Q-a-l/b+(N-l)b/2) 



Ci 



N 



here a = 01 + 02 + 03. Using the results of paper [67] where the correspondence between the 
integral J m (A, B, C\x) and differential equation of order m + 1 was established it is easy to derive 
the differential equations of order N — I + f for our correlation functions. 

4. The degenerate fields play even more important role if we consider our theory in the region b 2 < 0. 
For |6| < 1 parameter a\ in eq (14.51) becomes negative and we are in the other regime of the 
three exponential model. It is convenient to take — b 2 = -7- then as it follows from the results 
of [39j[63] for integer values of p > 2, three field model after quantum group restriction describe 
rational unitarian CFT M(N,p) with central charge 



"** = -WT7 [wrTi)) (5 ' 5) 

which is SU s^N+p-2)' ^ cose ^ CFT0 The conformal dimensions of the basic fields <\) k nm in this 
model are characterized by three integers k < N,m < p and n < p + N, with k = n — m (mod N) 
and are 

k(N-k) (pn-(p + N)m) 2 
m ' n ~ 2N(N + 2) + ApN(p + N) 1 j 

The fields (p^m U P ^° normalization calculated in [62] factor can be represented as <p\ m — > 
<Jke~ 2tuJmnip with 2u mn = b(n — 1) — (m — l)/b. The screening charges that are used to calcu- 
late the correlation functions in the M(N,p) theory have a form [68] 



It means that correlation functions in this theory can be derived from the correlation functions 
of the degenerate fields in parafermionic LFT by substitution b — > ib. In particular the structure 
constants in the M(n,p) theory can be expressed in terms of integrals calculated in Section 3. 

It was shown in ref [69] that the CFT models M(N,p) describe the critical behavior of generalized 
RSOS statistical models. The spin variable s in these models takes the integer values with the 
restriction \s — s'\ < N for the neighboring sites of lattice. These models can be considered at the 
random (dynamical) lattice. It is natural to think that the critical behavior of these "random" 
models is described by CFT M(N,p) coupled with parafermionic LFT. Solving parafermionic LFT 
on a sphere we will be able to find KPZ critical exponents and calculate the correlation functions 
in the coupled theory. To derive more interesting quantities as the partition function on the 
fluctuating disk we should solve boundary parafermionic LFT. One of the steps in this direction 
is the calculation of bulk-boundary parafermionic Selberg integrals which generalize the integrals 
calculated in [70] for the Liouville boundary CFT. 



13 



For non-integer but rational values of parameter p the model is rational but non-unitarian. 
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Appendix A. Parafermionic CFT (review) 

In this section we collect all needed facts about parafermionic CFT [20]. In particular, we give an 
evidence of the statement that correlation function in parafermionic CFT containing at most three 
order fields and arbitrary number of parafermionic currents is a polynomial (see (I2.19P ) and derive the 
set of properties ( 13. 5ft defining this polynomial uniquely up to a factor. We also explain our choice 
of normalization in ( 13 ,4p . In the end of the appendix we collect some formulae for the parafermionic 
polynomial ( I2.19P as well as for the integral ( 13 .ip in the dual case when parameters k±, k% and k% are 
related by the condition (12.2bl) . 

The Zn parafermionic algebra is generated by the holomorphic currents ^ k (z), k = 1, . . . , N — 1 of 
fractional spins 

a k(N-k) 



N ' 

The OPE of these currents can be written in a form 



A fe = 0. (a.i; 



2kk 

y k {z)y k ,{0) = C k>k ,z-^{y k+k ,{0) + O{z)), if k + k' <N, 

2(N-k)(N-k') 

#fc(s)tf fc /(0) = CV-^-fc' z w (^+^(0) + O(z)) , if k + k > N, (a.2) 
**(*)**-*(0) = (l + T(0) + 0(z*] 



where 



C ktk ' 



r(l + k + k')T(l + N- k)T(l + N-k') 



(A.3) 



_r(l + fc)r(l + k')T(l + N-k- k')T(l + N)_ 

Third line in ( ]A.2j) defines stress-energy tensor T(z) which generates the conformal algebra of the theory 
with the central charge 

2(JV-1) 



(N + 2) 

There is a set of antiholomorphic currents ^f k (z) with 



(A.4) 



A fe = 0, A k = k SE_A. (A.5) 

forming operator algebra similar to (IA.2p . Together they form the whole symmetry algebra of 
parafermionic theory and all other fields belong to some representations of this algebra. We note that 
the algebra (1A.2|) is generated by the first parafermionic current ^\{z) which we denote for simplicity 

as 

tf(z) = *i(z). (A.6) 
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Analogously we define ^f(z). It can be shown [20] that any field in theory can be obtained in the 
OPE's of the first parafermionic currents ^(z) and ^f(z) with one of the primary fields (it is convenient 
to include cr (z, z) = <tn(z,z) = 1.) 

<T k (z,z), k = 0,...,N-l. (A.7) 

Fields (Tk(z, z) are called by analogy with Ising model (which coincides with Z 2 theory) order fields and 
have conformal dimensions 

4 - 4 - 2N{N + 2 y (A.8) 

We note that flA.ll) . flA.2j) and flA.8j) have obvious symmetry with respect to k — > N — k. Sometimes it 
is convenient to define 

tt+CzO^ttjv-fcCz) and a+(z,z) = a N „ k (z,z). (A.9) 
Parafermionic theories enjoy Kramers- Wannier duality, so it is also reasonable to define disorder fields 

fi k {z,z), n+(z,z)= fi N _ k {z,z), k = 0, . . . , N — 1, (A. 10) 

which have exactly the same conformal dimension as order fields (1A.8D . It is convenient also to classify 
fields in theory with respect to both holomorphic and antiholomorphic Virasoro algebras^) 

$f M - ); q,q = k mod2Z, (All) 

where the parameters q and q are defined modulo 2N and 

(9.9) ^(q+N,q+ N) ^(q-N,q-N)- ^ ' 

The field $^ fcfc ) represents the order field a k , field ^_ kk \ corresponds to the disorder field fi k , 3>( 2 fco) 
corresponds to the parafermionic current ^ k while ^ 02k ) corresponds to ^/ k (we identify ^0 = ^jv 
with identity operator). The system of the fields ^ gk ) obtained by the application of the holomorphic 
part of parafermionic algebra to one of the order fields can be conveniently drawn as shown on fig 
[TJ Parafermionic currents ty k are drawn by big dots, order fields a k by circles and disorder fields \i k 
by crosses. All other fields $( k+2 i k) w hh I = {—2k, . . . , N — k — 1) shown by small dots on fig [T] are 
parafermionic descendants of a k or fj, k with conformal dimensions 

A($f fc+2W ) = 4 + KN ~* " , A($f fc _ 2W ) = d k + (A.13) 

where d k is conformal dimension (|A.8[) of the order field a k . We note that parafermionic current ^(z) 
acts horizontally from the left to the right on fig [TJ The fields local with respect to themselves and with 
order fields are Jjj) 

One can consider the multipoint correlation function 

& (n\ fci,fa,fc 8 )( a . lj X2i X3 \t u ... : t n ) = (a kl (x u x x ) <T k2 (x 2 , x 2 )cr k3 (x 3 , z 3 )tf (ti) ■ ■ ■ #(*„)#(*„)), (A.14) 



14 We note that the field primary with respect to Virasoro algebra is not necessary primary with respect to parafermionic 
algebra (|A.2|) . 

15 Analogously the fields _ q s are local with themselves and disorder fields fi k - 
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Figure 1: Digram describing the system of fields *&( q fy m the Z w model. 



where n = mN — k. It is clear that the sufficient condition that this function is non-zero follows from 
the similar condition for the four-point function 







i Nlkl ' k2 ' k3 \xuX 2 ,x 3 \t) = {a kl (x lj Xi)a k2 (x 2j x 2 )a k;i (x 3 ,x 3 )^ N ^ k (t)^ N _ k (t)}. (A. 15) 



In order to find out whether correlation function ( 1A.15I) is non-zero we consider the operator product 
expansion [20] 



E 



s 

1=0 



(91,91) (92,92) 



(jjfcl+&2— 2i 

(91+92,91+92) 



E 



1=0 



tf ) k 1 +k 2 -N+2l 

(91+92-iV, 91+92-A^) 



s = min(A; 1 , k 2 ) for ki + k 2 < N, 
, s = min(N - k u N - k 2 ) for k Y + k 2 > N. 



(A.16) 

Evidently ( fA~T5l) is non-zero if and only if in the OPE a\a 2 a 3 there is a field ^> k ^k- Analyzing ( 1A.16I) 
one can arrive that this is possible if one of the following conditions 



kx + k 2 + k 3 = 2k, < kj < k < N, 

h + k 2 + k 3 = N + 2k, < k < kj < N, 



(A. 17a) 
(A.17b) 



is satisfied. We consider in details the condition (IA.17aj) . In this case correlation function (IA.14[) can 
be written in a form 



~(h h h \ \-pi N \ k i>k2,k 3 ) 



j 

I X 12 1 25 i2 I Xl3 1 2<5i3 | X23 1 25 23 ^ |£. _ rj.^. |^ 



N 



(A.18) 



where a:^ = Xj — Xj, = — t,- and 



fcjfcj (iV - fc + l)(fci + fcj - k) 
ij ~ N{N + 2) 



(A^ + 2) 
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The normalization factor H(&i, fa, fa) in (1A.18j) will be denned below. Clearly P(xi\tj) denned by 
(1A.18P is a polynomial symmetric in ti,...,t n . Projective invariance of correlation function (1A. 18|) 
requires that this polynomial satisfy the set of properties (j3.5aj) - (l3.5cp . The last two properties in 
( 13. 5 p follow immediately from the operator product expansions (IA.2[) and (jA.16|) . Namely, correlation 
function (1A. 14|) should respect OPE's 



N N-kj 



(A.19) 



From other side applying (IA. 19[) to (IA. 18|) one arrives to unwanted poles in (|A.18[) . The condition that 
these poles are cancelled is equivalent tco 



n-N 



X 



N 



n+kj-N 



(A.20) 



P£ lkl ' k2 ' k3) {x 1 ,x 2 ,x 3 \x :j ,...,x j ,t 1 ,...,t n+kj - N ) ~ \ \ (ti-xj), for j = 1,2,3, 



N-k, 



i=l 



Clearly eq (1A.20I) implies (13.5d|) and fl3.5e[) . One can see that there are no more obstructions and 
correlation function (IA. 18[) with appropriate polynomial satisfy all eligibility requirements. Moreover, 
as it will be proven in appendix [B] the set of properties (13. 5ft determines the polynomial P(xi\tj) in 
(IA.18P uniquely up to factor which can be defined as follows. In order to obtain the four-point function 
(1A.15P from the multipoint one (1A.14D one has to fuse n = mN — k points t\, . . . ,t n into (m — 1) 
clusters with N points and one cluster with N — k points taking into account the structure constants 
of parafermionic algebra (IA.3j) . i.e. 



N N N-k 





m—l 


~N\(N-k)\~ 






N N ~ k k\ 



N-k- 





in ] 



In our conventions we defined normalization of the polynomial P^' kl ' k2 ' ka \xi, X2, £3^1, 
according to dA~20D and ([Oil as 



(A.21) 



t n ) in ([31D 



p(AT|fel,fc 2 ,fe3)/ 





m—l 


~N\(N-k)\~ 


IN*. 




N N-k k ] 



X3\ yi,---,yi , y m -i,...,y m -i , t^^t) = — - , , .... x 

N N N-k 

(m—l) m-1 

x II fa-fc) 2 " n(yj-*i) h (yj-*J*(yi-z*) h '(yj-t) 3{N - h) - ( A - 22 ) 



i<j=i 



3=1 



According to (IA.2ip and (1A.221) the normalization constant S(fei, k2, fa) in (IA. 18[) is defined through 
the four-point correlation function (IA. 15[) as 



0it 



(N\k 1 ,k 2 ,k 3 ) (rr u\ 
k [Xl,X2,X 3 \t) 



E(h, fa, fa) 



\X12r 012 F13 



2<5i 



F23 



2<?23 



2(N-h)kj ' 



(A.23) 



16 



Here sign ~ means proportional up to some polynomial. 
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From operator product expansion OA. 160 one immediately concludes thatl 17 ! 

E(kx,k 2 ,k 3 ) = K^^t-fc^-*!-^ = 

— /cb kl cb N ~ k2 cb N ~ k3 \/<b N ~ k3 (b k s (b N \ (A OA) 

" \ (fci, fci) (ife-iV, k 2 -N) ^(N-ki-ka, N-ki-k 2 )' \ (fci+fe-Af, fci+fc 2 -iV) (fc3, fc 3 )^(iV-2A:, N-2k)h j 

An explicit expression for the three-point functions in Zjy parafermionic theory was obtained in 



2ii $ 2 J2 _ $ 2j 3 



(2mi,2mi) (2rri2,2m2) (2m3,2m3) 

where first two factors are 3j-symbols and 



> 



Jl 32 33 

mi m 2 m 3 



3i 32 33 

777.1 m 2 "73 



(A.25) 



t(j'l,J 2 , is) 



7 



1 



N + 2 



7 



N - 2ji + 1 



7 



jV - 2j 2 + 1 
N + 2 



7 



AT + 2 



x 



n(ji + 32 + 33 + i)n(j 2 + j 3 - jjMh + j 3 - jzMh + j 2 - j 3 ) 

n(2 jl )n(2j 2 )n(2j 3 ) 



(A.26) 



with 



J 



N + 2 



n(ife) = fc! JJ 7 

3=1 

In (1A.25|) sign ~ means that the l.h.s. is not necessary zero if the r.h.s. does. One has to use reflection 
properties ( 1A.12|) for the fields 2 m) * n AA.25I) to ensure that the r.h.s. is non-zero (if it can not be 
done then this correlation function is indeed vanish). In our case we tuned parameters in (1A.24|) using 
this prescription and get 

kUN-k\\ 

(A.27) 



E(k u k 2 , h) = k ' (yN AT . ^' p(h, k 2 , k 3 ), 



where p{k\, k 2 , k 3 ) is given by (12. 4p . With explanations above we think that the choice of normalization 



is 



in (13 .4p becomes clea 

In the end of this appendix we give an explicit expression for the parafermionic polynomial (I2.19p 
in the dual case when parameters fci, k 2 and k 3 are related by 

h + k 2 + k 3 = N + 2k, < k < kj < N. (A.28) 

In this case polynomial defined by (I2.19P is expressed in terms of polynomial Pm kl ' h2 (ti, ■ ■ ■ ,t n ) as 
(here again n = mN — k) 

^ l * 1,fca '* s) (ti,...,t„) = (-l)^^ 1 ) 



k\{N-k)\ / 
— p(k 1 ,k 2 ,k 3 ) 



Pf" fe ' y (ti,.,g, (A.29) 



where p(ki, k 2 , k 3 ) is given by (J23D and P^ kl ' ha,ka) ( i t 1 , . . . ,t n ) is defined similar to ([23]) by 
Pf 1 **'^!, ■■■,t n ) = {-l) m ^~^ lim x\~ mk * PW^MM)^ hx3 \ tu . . . )tn ). 



£3— >oo 



(A.30) 



17 In the second line in (IA.24|) we used the relations (|A.12[) . 
One has to take a simultaneous look at eqs (fQ8|) . (gUgi and (gU) to realize our choice of normalization. 
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Polynomial Pm kl ' k2 ' k3 \x\, x 2 , x 3 \ti, . . . ,t n ) satisfies the set of properties similar to (13. 5p 

P(Xxi\Xtj) = X m ' N ~ k P(xi\tj), (A.31a) 



3 n 



3=1 q=l 

P(xi + X\tj + A) = P(xi\tj), (A.31c) 

P(X 1 ,X 2 ,X 3 \X, . . .X,t U . . -tn-N-l) = 0, (A.31d) 
N+l 

P(xi,X2,x 3 \x j ,...x j ,t 1 ,...tn+k j -N-i) = for j = 1,2,3. (A.31e) 

N+l-kj 

One can shown that (1A.31|) determine the polynomial Pm^ kl ' k2 ' k3 \xi,x 2 ,x 3 \ti, . . . ,t n ) unambiguously 
up to a total factor. It can be written similar to (13. 6 p as 

P? kM \x u x 2 , x 3 \t u ...,t n )= A" 1 S\p( Xl , x 2 , x 3 \h, t n )}, (A.32) 

where A m is given by (13. 7p and the basic "monomial" p(xi,x 2 ,x 3 \ti, . . . ,t n ) is defined as follows. We 
divide our points ti,...,t n into N groups 

(ti,t 2 , ■ ■ ■ , t n ) — > (si, . . . , s^), 

where N — k groups contain m different elements and the remaining k groups contain m — 1 different 
elements. All groups (s±, . . . , sn) are divided into four sets with ki — k, k 2 — k and k 3 — k groups (m 
elements in each group) and k groups ((m — 1) elements in each group) 

Xl — ( s l> • • • j Ski-k), X2 = (s kl -k+l, • • • , SN-k-s), 

X3 — (sjv-it 3 +i) • • • , $N-k), Xi — (SjV-fe+1) • • • ) $n)- 
Then the basic monomial is 

p{xi,x 2 ,x 3 \ti, ...,£„) = 

N 

= n n & - i tf n - x ^ n & - ^ n & - ^ n - x ^ - x ^ - ^ ( a - 33 ) 

q=i i<j£ Sq jexi j'gx2 iex3 jex4 

We can define also the dual integral 

i i 

„ „ n n 

j(N\kuk.M) {ABlg) = j jY[ t A {1 _ t . )B "Q \ ti - tj \^P^ k ^\t 1 ,...,t n )dt 1 ...dt n , (A.34) 







i=l i<j=l 



where the integers ki, k 2 , k 3 and k are restricted by the condition (IA.28I) . It can be calculated using 
the method described in section [3j The result is 

l l G^! ks+k (l + A + B + ((m + P )N-k-l) g + m + p) 



gKV + A + (m - ljjVg) G[ N 2 : k \m + B + (m- l)Ng) 
G { "-£(2m + A + B + ((2m - 1)N - k - l)g) 

where function G k N \x) is defined by (13.161) and function Cm' k \g) is given by f!3.35j) . 
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Appendix B. Proof of proposition 



EL 



In this appendix we give a proof of the following statement used in section [3j Namely, the polynomial 
P(xi\tj) defined by the set of properties (13.51) exists and unique up to a total factor. The existence was 
demonstrated in section [3] by explicit expression, while the uniqueness will be proven below (we follow 
mainly [54"]f*^l. 

Let us consider the space W of symmetric polynomials P(ti, . . . , t n ) satisfying the following condition 

P( x, . , x , ti, . . . , t n _jy_i) = 0. (B.l) 

N+l 



Proposition B.l For any polynomial P(ti, . . . , t n ) E W 



P[Xj . . . , X, t\j • • • , tn—N 

)~l[(t J -x)\ (B.2) 

n i= l 



where sign ~ means proportional up to some polynomial. 

Proof @ Let us take in P(ui, . . . , u N+ i,t 2 , . . . , i n -Jv) 



Uj = u + A£j, with u = 1 N+1 and £j = 0. (B.3) 

+ 3=1 

We can expand at A — > 

P(w + ACi, . . . , w + A^iv+i, *2, • • • , t n . N ) = P {0) (u\t 2 , . . . , t n ) + A Sl (0P (1) (M|t 2 , . . . , t n ) + . . . , (B.4) 

where s 3 -(£) are the elementary symmetric polynomials in (^ 1; . . . Due to the conditions (IB . 1 j) 

and flB.3j) first two terms in the expansion (IB.4[) are actually vanish. So that the expansion (1B.4|) takes 
the form 

P(u + ACi, ...,« + A^iv+i, * 2 , • • • , in-jv) = A 2 s 2 (0P (2) Ht 2 , • • • , t n ) + 0(A 3 ). (B.5) 

Now we take u\ = ■ ■ ■ = un = x and u^ + \ = t\ which is equivalent to u — (Nx + ti)/(N + 1) 
& = ■ ■ ■ = £ N = ( x - h)/(X(N + 1)) and = ~N(x - *i)/(A(iV + 1)) and find that 

P(x, ,x,t u ..., t„_Ar) ~ (ti - x) 2 , 

N 

which due to the symmetry reasons implies (1B.2j) . 

19 Along the same lines one can prove uniqueness of the dual polynomial defined by the set of properties (|A.3f j) . 
20 We note that there is a simple explanation of proposition IB . II from different point of view. Namely, one can consider 
dual space W* which is the space of all n— point correlation functions 

P(tx,...,t n ) = {i)\e(h)...e(t n )\0), 

of the operator e(t) = J>2j ejt - -* -1 , [e<,ej] = satisfying e(t) N+1 = 0, where runs over W*. Clearly the relation 
e(t) N+1 = implies e'{t)e(t) N = 0. It follows from the relation e{t) N+1 = that P{x, ...,x,h,.. . ,t n - N -i) = 0, while 

N+l 



N+l 



from e'(t)e(t) N — we find that P(x, . . . , x, tx, . . . , tj, . . . , t n _jv) 

JV 



= 0, which is equivalent to (|B.2[) . 
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Let us define the space Uf labeled by Young tableau Y with n cells 

y = (pi s P2>---)> $^^ = n ' 

as the space of polynomials satisfying 

P eU 9 if P{Y') = for any Y' > Y, (B.6) 
such that YuPj = n - I n AB.6I) the notation P(F) means 



P(Y) = P(|^_^, - - - ), (B.7) 

Pi P2 

and by definition partition Y' is greater than partition Y if there exists j such that p\ = pi for % < j 
and p'j > pj. This ordering is called lexicographic. It is clear that W = U-p where Yq is the following 
partition (pi = N + 1, p% — 1, . . . ) 



N + l 



(B.8) 



n-N 



It is also clear that Uf C C/yf, if K < V . Thus we have a filtration 

C • • • C U 2 C C/i C C/ , 



(B.9) 



where we denote U$ = W, U\ = U-g- with Y\ being the partition with n cells preceding the master 
partition Y in lexicographic order and so on. Let us define the quotient space 



U q = U q -!\U q , 
which can be identified with the space of polynomials 

P(Y q ) ^ and P(Y) = for any Y > Y q . 



(B.10) 



(B.ll) 



As for any polynomial P G U q P(Y q ) / we can restrict it on the diagram Y q = (p[ ,p 2 > • • • ) an d 
define according to <\B.7\i the polynomial XqiVii V2, ■ ■ ■ ) as 



Xq(Vl,V2,>") = P{Y q ). 



(B.12) 
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It is clear that there is one to one correspondence between polynomials P G U q and polynomials 
x{yi:V2, • • • ) defined by (IB. 12j) . i. e. the map (IB.12P is injective. The polynomial 2/2, • • • ) should 
vanish at the diagonal i/i = yj as 



XqiVi, V2, ■ ■ ■ ) ~ {Vi ~ yj) 2 ^ for i < j. (B.13) 



Indeed, for any polynomial P G U q we can consider polynomial g(ti, . . . ,t 



g(ti, . . . ,t M (,)) = P{. . .,yi-i, ■ ■ ■ ,y i -i,t l , . . . ,t ( q ),y i+1 , . . . ,y i+ i, . . ., 

Pi +Pj v w ' Pi v w ' 

„(«) „(«) 

Pi-l Pi+l 



Vj-ii • • • > 2/j— is > • • • > Z/j+i) • • • j Z/i+i) • • • ), ( B -14) 

V / Pi ^ 1 Pi ^Pj v „ / 



„(<?) (8) 



i.e. we specify polynomial P on almost whole diagram Y q except i-th and j-th columns (here z < j and 
hence p- 9 '' > pf )■ Clearly, the polynomial g(ti, . . . ,t <,) (,)) defined by (IB. 141) is non-zero. Moreover, 



it should vanish if p[ + 1 variables t{ coincide, because in this case we get polynomial P specified on a 



diagram which is definitely greater than Y q . Now as follows from proposition lB.il 

Pf 

g(yi, ■ • • , Vi, h, . . . , t w) ~ Ylivi - t a ) 



\2 



Pi 



(«) a=1 



and hence 



eiVi, ...,yi,yj,..., yj) ~ {yi - yj) 2p{ * } - (B.15) 



p^ p$ 9) 



It follows from (1B.15|) that 



x q (yuy2,---)-U^~y^ q) (B.16) 

Now we restrict the space W and define the space 

yyikitete) c W as the space of polynomials 

P G W satisfying 

) = for 0<A; a <iV a = 1,2, 3. (B.17) 

N+l-k a 

Where points x\, £2 and X3 are some special points in a sense that we allow the coefficients of the 
polynomial P depend on x a . In a same way we define the quotient spaces 

Proposition B.2 For any polynomial P(t u . . . , t n ) G W^ klM ^ 3) 

3 

P(x, . . . , x, h, . . . , t n „ N+l ) ~ - x a )^, (B.18) 

iV-Z a=1 



21 



It depends also on variables {y\, . . . , . . . , 2/j-i, Vj+i, ■ ■ ■ ), but this dependence is hidden for shortness. 
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where 



k a I if / < k a , 
if I > k a . 



Proof Let us denote for simplicity f(ti, . . . , ijv-z) = P{ti, ■ ■ ■ , fjv— i, t^-i+i, ■ ■ ■)■ Using the expansion 
valid for any symmetric function f(t\,... , tjv-z) 







f(x + • • • , x + \£ N -i) = f{x, ...,x) + \^2 £j-^-f(h, x,...,x) 



d 2 J d 2 



+ 



tl =x 



f(ti,h,x, 



,x 



tl=t2=X 



+0(A 3 ), 



and taking x = x a we immediately conclude that due to (IB. 171) 

P(x a + \£x, ■ ■ ■ , x a + X^n-i, h,..., t n - N+ i) ~ \ Ka 
taking ^ — (x — x a )/X we arrive at (1B.18j) . 

It follows from proposition IB. 21 that the polynomial Xq{Vii Vi-i ■ ■ ■ ) specifies even more 



(B.19) 



,2p 



(«) 



i a=l 



i<j 



where 



K 



(?) 
j,a 



pf +k a -N if V f > N — k a 
if p$ 9) <N-k a . 



(?) 



Now let us define according to (1B.20I) 



^-0 + 4>( 2 £p. w + X>S) 

Pj \ i>j a=l J 



(B.20) 



(B.2i; 



which is the minimal possible degree of polynomial XqiVii 2/2; ■ ■ • ) i n variable yj divided by (we note 

is at least 

(B.22) 



that vf^ is not necessary an integer number). It is clear that any polynomial P G U q 
of degree 

u {q) = [maxz/^1 

in each variable t\, . . . , t n where \x~\ means ceiling function i.e. \x] = min{m G Z|m > x}. 

Now we turn to the special case of polynomials defined by the set of properties (13. 5p . In this case 
integers k\, k 2 and k% are restricted by the condition (I2.2ap and n = mN — k. It follows from (13.5bj) 
the polynomial P^ kl ' k2 ' k3 \xi, x 2 , x$\ti, . . . ,t n ) is of degree 2(m — 1) in each variable tj. We will show 
below that 2(m — 1) is the minimal possible degree for the polynomial from W^ kl,k2 ' k ^ with the condition 
that integers k, ki, k 2 and k 3 are related by (12.2aj) and n = mN — k and if such a polynomial exists 
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it should be unique. Let us compute Vj denned by (IB. 211) . By definition the diagram Y\ preceding i 
lexicographic ordering the master diagram Yq drawn in ( IB. 81) is 



N 



••• ••• ••■ 



1111) 

• • • • • • 1 



m 



(N-k) 

I 



(B.23) 



Using definition (IB.21[) we find 




Now, using the condition 



D + 



D + 



N ' 
(i) 



for j 



1,.. 



N-k ' 



m. 



5> 



2k, 



for j 

0<k a <k<N, 



we find that K r 



and hence 



Moreover, one can easily see tha 



(i) 



v 



(i) 



V, 



(1) 



2(m - 1) 



, m 



[B.24) 



(B.25) 
(B.26) 



„(«) > U W f or q > i. 

So, if there exists a polynomial P G jy( fcl ' fc2 > fc 3) with ^ k a = 2k and < fc a < < of degree 2(m — 1) 
in each variable t±, . . . , t n it should be unique, because in this case dim(f/ 1 (fcl ' fc2 ' /C3) ) < 1, while as follows 
from flB~26l) dim(U^ klMM) ) = for q > 1 and hence 



dim(W (kl 



^dim(^ fcl - 



< 1. 



(B.27) 



An explicit form of this polynomial was given in section [3j Of course, the total normalization remains 
unfixed under above consideration. We chose it regarding our needs as in (13. 6p . 



22 It is clear that for any diagram Y q < Y\ with the total number of columns greater than m the degree v^> > 2(rn — 1) 
(it follows immediately from the definition (IB.21[1 ). So, we must consider only diagrams with m columns. Let us take one 
of them then evidently p^-i — N — I and = N — k + I + e with some positive integer I and non- negative integer e. 



From definition (|B.21I) we find that 



,(9) 



2(m - 1) 



(N-l) 



4 ;-2fc + 2e + ]>>^_ liQ 



Using that J2l=i K m-i a — El=i i^a — = 2fc — 3/ we find that the expression in brackets is greater than I + 2e and hence 
>2(m-l). 
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Appendix C. Proof of proposition 13. 2L 

In this appendix we give a proof of proposition 13.21 For simplicity, we consider the case k\ — h% = 
k 3 = 0. We define the polynomial P(t\, . . . ,t n ) with n = mN which is proportional up to a factor 
to the polynomial P^'°'° (ii, . . . ,t n ). It can be represented as a symmetrization of the "monomial" 
p(ti, ...,t n ) 

N 

P (t l ,...,t n )=n n (u-trf, (c.i) 

q=l i<jes q 

where (si, . . . ,sn) are the groups of m different variables tj. For non-negative integer values of the 
parameter g we can define the polynomial 

£>(*!, ...,t n ) = Y[\U- t^Pih, t n ), (C.2) 

i<j 

which can be expanded into the sum 

D(t U ...,t n ) = J2 C ^« t l 1 --- t n- ( C - 3 ) 

The proposition 13.21 in this case reads as 

Proposition C.I For V\ < v 2 < ■ ■ ■ < u n the non-zero values of c Ulj ,,, tUn in ( 1C.3I) occur when the integers 
z/ p AT +J with p — 0, . . . , m — 1 and j — 1, . . . , N satisfy the inequalities 

p+(pN + j- l)g < is pN+j <(m + p-l) + ((m + p)N +j- 2)g. (C.4) 

Proof Clearly, the polynomial (1C.3P is a homogeneous symmetric polynomial of degree 

(m — l)n + n(n — l)g, 

i.e. X]j=i = ( m — l) n + n ( n ~ l)fl l - This summation formula together with the assumed ordering 
vi < ^2 < • ■ ■ < implies that 

i/ n >m-l + (n-l)p. (C.5) 
To estimate z/ n _fc for we use the fact that 

n—k 

> (m- l)n + n(n- l)g - z/f", (C.6) 

j=i 

where z/™ ax is the maximal collective degree of the polynomial D(ti, . . . , t n ) in variables t n -k+i, ■ ■ ■ ,t n . 
In order to find z/™ ax we take 

t-n-k+l = r £n-fc+l) tn-k+2 = T £,n-k+2, • • • = r £nj 

and find the leading asymptotic of D(t±, . . . , t n -k, r^ n _jt + i, . . . , r£ n ) in the limit r — >■ oo. First of all, we 
note that 

max -max i " max 

V k — V k k ' 
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where the degree z/™ ax comes for the asymptotic of the polynomial Y\.i<j — \ 2g wnne ^™ ax comes from 
the polynomial P(t\, . . . ,t n ). The first degree is a simple one 

z>, max = k(2n - k - l)g. (C.7) 

In order to find z>jf ax we note that the asymptotic of the polynomial P(ii, . . . , £ n -fc> r £n.-fc+i) • • • , r £n) in 
the limit r — > oo depends essentially on an integer part of the fraction k/N. It is convenient to represent 

k = (m — p)N — j, p = 0,...,m-l, j = l,...,JV, (C.8) 

i.e. n — k = pN + j (we remind that n = mN). It is clear that the leading asymptotic of the polynomial 
P(t\, . . . , t n -k, r£ n -k+i, • • • , T^n) comes from the basic "monomials" such that j groups contain (m— p— 1) 
variables r^j and (N — j) groups contain (m — p) variables. Using simple algebra one can find 

v™ ax = (m-p)(m + p-l)N -2pj. (C.9) 

Using (\QM), (jCJ]h (1081) and ([09]) one can find 

pN+j 
i=l 

and hence from the assumed ordering of the integers V\ . . . v P N+j we get 

p(N - j) 



VpN+j >P + (pN + j — l)g 



pN + j 



We note that the number p(N — j)/(pN + j) is always smaller than 1 and hence it can be dropped in 
the inequality above. So, we get a final formula 

"pN+j >P + (pN + j- l)g, 
which is the left inequality in (]C.4j) . Now, using the fact that 

n 

D(h, ...,*„) = Htf m - 1)+2in - 1)s D(l/ tl , l/Q, (CIO) 

we can expand the r.h.s. of ( 1C.10I) as in ( 1C.3I) with v[ < v' 2 < • ■ ■ < v' n . We note that ( 1C.10I) implies 
that v[ = 2(m — 1) + 2{n — l)g — is n +i-i and hence the inequality v' P N+j > P + (pN + j — l)g implies that 

v pN+j <(m + p-l) + ((m + p)N + j- 2)g, 

which is the right inequality in ( 1C.4I) . 

In the case of general parameters kj and k the proof of proposition 13.21 can be performed in a same 
way. We hope that the explicit form of the inequalities (13 . 13[) will help interested reader to repeat above 
arguments in this general case too. 
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Appendix D. Sigma model representation and minisuperspace 

approximation for three exponential model 

In this appendix we consider minisuperspace approximation for CFT (14. 3[) and its sigma model repre- 
sentation. For this purpose we consider this theory in the dual representation with screening fields Jj 
defined by (I4.27P and ( I4.28p . At first we consider this model with only one screening field J±Ji, which 
corresponds to the parafermionic Liouville CFT. It follows from the explicit form of the current J\ that 
the action 

A 1 = J ^f^B^-^M^+^j (D.l) 

where our dilaton field is $ = + ^"V^i + 20^2 can ^ e wr hten in terms of sigma model with 

complex metric. In minisuperspace approximation when iV and parameters a, oti are large we can use 
the relation: a 2 = a\ + a 2 . To simplify the notations we denote <pi = Xi,ip 2 = x 2 ,(p 3 = x o- Then our 
metric can be written in the form. 

x 2 x 2 x 2 

ds 2 = (1 + ma^^jdig + (1 — ma|e^ )dx\ + (1 — ma^ )dx\— 

x 2 x 2 x 2 

— 2iaaime 2a 2 dx^dxi + 2iaa 2 me 2a 2 dx^dx 2 + 2a\a 2 me 2a i dx\dx 2 (D.2) 

where parameter m = 47r/z. Our metric qq. is complex in the "Euclidean" space, however in "Minkowski" 
space: xq — > ixq this metric is real regular metric with signature (—1, 1, 1). This metric possesses the 
properties g = det[gik] = l,g ik (m) = gik(-m). The metric is not conformally flat. It has non-zero 
Ricci tensor. The scalar curvature for this metric R = — ^me 2a 2 . It is small for x 2 — > —00 and becomes 
large for x 2 — * 00, where we expect strong interaction. To compensate the Ricci tensor we should take 
into account dilaton field $ = j^Xq + ^-^i + 2 ^ x 2- ^ ne ac tion for renormalization group flow is 



S = J d z Xyfge-*{R+\V®\ 



and the condition of conformal invariance has a form 

R ik + V*V fe $ = (D.3) 

It can be checked that this equation is satisfied with our linear dilaton field $. The equation that 
describes zero modes dynamics in minisuperspace approximation has a form 

- ^=d ie -*^gg ik d k Q = A0 (D.4) 

We see that our metric depends only on the variable x 2 , it follows from the fact that field J\J\ interacts 
essentially only with field ip 2 . It means that function depends non-trivially only on variable x 2 . It is 
natural take this function in the form 

= exp{$/2 + iqxo + ipiXi}^/(x 2 ). 



Then eq ( ID.4j) can be reduced to the hypergeometric form. The unique regular at x 2 — > 00 properly 



normalized solution to this equation with eigenvalue 



16a 2 16af Ifta^ 
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has a form 



r(| -2ia 2 p 2 + 2aq-2ia 1 p 1 )T( 1 - 2ia 2 p 2 - 2aq + 2ia 1 p 1 ) 2 ,-. 2iaoVo +i 

W = — — - (man) zt " 2 ^ 2+ 2 

Y{-Aia 2 p 2 ) V 2J 



x 2 iatpx <*q\ { 1 
— + x 2 { U F [A, B, 1, j 



x exp <| + x 2 ( — ^ "Z-^ e "^ ) ( D - 6 ) 



where A = ^ + 2ia 2 p 2 + 2ag — 2z«ipi, £> = ~ — 2ia 2 p 2 + 2aq — 2ia\p\ and numerical factor in eq ( ID. 61) is 
chosen in such way that asymptotic of function \l/(x2) in the free region x 2 — > oo has a canonical form 

*(x 2 ) = e^ 2 + 5 2 (p 1 ,p 2 , g)e"^' 2 (D.7) 
It follows from explicit form of our solution that 

5 2 (pi,P2,g) = (ma 2 ) 2P2 T^^ v /i , r ^ V ^ D - 8 ) 

1 {-4:ta 2 p 2 )'y{^ + 2ia 2 p 2 + 2aq - 2iaipi) 

If we take into account that m = 47r/z and that as follows from eq (14.321) 

S)K W (35|)) 2 ^ (4 ^ (D - 9) 

we see from eqs ( I4.19|l4.2ip that in the minisuperspace limit where a,i » 1 and 

1 1 

2ajdj = - + i2ajPj, 2cra 3 = - + 2aq 

are finite S 2 (pi,p 2 ,q) coincides with reflection amplitude R 2 {a). 

In the same way if we substitute in the action (ID.lj) the term /iJx Ji by the term J1J 2 J 2 and in the 
eq f ]D.9j) ji — >■ /i and a 2 — >■ «i we can calculate coefficient Si(pi,p 2 ,q) for outcoming wave e~ ipiXl and 
check that it coincides with minisuperspace limit of reflection amplitude Ri(a). 

To derive the minisuperspace limit of the reflection amplitude i? 12 (a) we should take in the action 
both terms [iJ\J\ + J1J 2 J 2 . However, in this way we have some problem with the correct definition of 
the action. The reason is that the operator product expansion of currents J\ and J 2 is rather singular 

Mz) M z') = (, - ;<)- 2 exp + ^ ) + (D.IO) 

and we should renormalize our action by adding counterterms in such a way that the symmetry of our 
action becomes not broken. For this we note that the primary field with dimension 

G(z) = (iad z ip 3 + a^ipi + a 2 d z Lp 2 ) exp | _ z ^ 3 ( z j (D.ll) 



2a 

commutes with all screening charges S+ + , Sl_ and Sl + of three field theory H4.3[) and plays the role of 
the generator of the symmetry of this theory. It eommutes also with screening operators Q, = fj,(r.)<iz. 

It is easy to check that field G(z) commutes also with screening operators Q\ = <j) Li(z)dz generated by 
the currents 

h(z) = —^==(iad z <p3 + atxd s <pi + a 2 d z ip 2 ) exp (D-12) 
V4iV V 2a 2 / 
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I 2 {z) = —^=(iad z ip 3 + a x d z ^ x + a 2 d z (p 2 ) exp ( fl^l ] (D.13) 
V4iV V 2a i / 

It is not surprising because currents Jj differ from currents J« by total derivatives which do not contribute 
to the screening operators. The main advantage of currents Ii is that their operator product expansion 
is not singular. That permits us to write well defined action in the form 

r ( 1 3 ~ - i? (2) \ 

A n = d2x \^T J2( d ^i) 2 ~ N ^hh - NJihh + -— $ (D.14) 

J \ 1D7T 07T J 

with the same parameters /i, fx and dilaton field The metric for sigma model corresponding to this 
action will have a form 

ds 2 = (1 + a 2 U)dx 2 + (1 - a$U)dx\ + (1 - a 2 U)dx\- 

— Haa-JJ dx§dx\ — 2iaa 2 U dxodx 2 — 2a\a 2 Udx\dx 2 (D.15) 

where = me 2a ^ + me 2 °i with m = 47r/x, m = 47r/x. 

It can be checked that this metric satisfies the conformal invariance condition 

^ + V,V fe $ = 0. (D.16) 

The curvature of this metric is 

R = __( me 2^" + me^T) = —U(x). (D.17) 

It means that interaction is strong for large positive X\ )2 . It is convenient to modify equation (ID.4[) for 
the wave function describing zero modes dynamics and to write it in the formal 



{-^y/gfdu + ^j^ = {p\+pl + g 2 )H. (D.18) 

The solution of this equation differs from the solution to eq (ID. 41) by the trivial factor exp[$/2]. The 
solution of our equation can be represented in the form 

Z(x ,x u x 2 ) =e iqx °-V(x l ,x 2 ) 
where function ^(xi,x 2 ) satisfies to the following equation in the partial derivatives (di = d xi ; d\ = d X2 ) 
_ dl _ U( x ) ^ aidl + a 2 d 2 ) 2 + ^(1 - 4ag)<9 1 + ^(1 - Aaq)d 2 )) *( Xl ,x 2 ) 
(tq '{Aaq - 1) + -M U(x) + {pl+p 2 2 ) j #(zi,rr 2 ). (D.19) 



2 v " 16, 
To solve our equation we substitute our function as 



ty( Xl ,x 2 ) = -j^- J dqidq 2 ^(q 1) q 2 ) exp (-iq lXl - iq 2 x 2 ) . (D.20) 



2:S 



The application of this equation to the analysis of the non-trivial Ricci flows can be found in [39 1I7TH73 
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Then after simple transformation eq (ID. 19|) can be written in the form 



V 2 



2iaiq± + 2ia 2 q 2 + 2aq + - m^(q 1 - - — , q 2 ) + m^/(q 1 , q 2 - 



2 J \ y ^ 2ai * ' ^ 2a 2 

= 4 ({<£ - p\) + (qj - pi)) *( qi , q 2 ). (D.21) 

The unique solution to this equation corresponding to the regular function ^/(xi,x 2 ) at Xi >2 — > oo is 

<f lf \ n , T{2ia 1 (q 1 -p 1 ))T(2ia 1 {q 1 +p 1 )) 

y{qi,q 2 ) = D(p u p 2 ,q) — — — — , x 

7(2«aigi + 2ia 2 q 2 + 2aq) 

T{2ia 2 {q 2 - p 2 ))T{2ia 2 {q 2 +p 2 ))(m?)- i2n « 1 (mQ^ ,2fl2 « 2 

where the constant 



D = (4a 1 a 2 )(ma 2 )" aaipi (m«2 



2 ^ i2aipi ( ^^-a a2P2 -y(-2iaipi - 2ia 2 p 2 + 2aq) 

T(-4iaip 1 )T(-'iia 2 p 2 ) 



The contours of integration in the eq (1D.20l) go along the real axis avoiding the poles at = ±pi 
counterclockwise. At X\^ 2 — > — oo (free region) we can close the contours in the upper half plane. The 
main term of the asymptotic of the wave function will be defined by the poles at q,i = ±pi and will have 
a form 

x 2 ) = e ipiXl+ip2X2 + S ie - ipiXl+ip2X2 + S 2 e ipiXl ~ ip2X2 + S lt2 e~ ipiXl - ip2X2 (D.22) 
where coefficients S 2 (pi,p 2 ,q) and Si{j>\,p 2 ,q) were calculated above and coefficient Si )2 is 

r(-4?aipi)r(-4«a 2 p2) 
7(-2i«ipi - 2ia 2 p 2 + 2aq) 
7(+2i«iPi + 2ia 2 p 2 + 2aq) ' 

If we take into account that m = AiiJl and that = f^TA* fi^J J — (^u4a^) 2ai we derive that 

Si 2 coincides with minisuperspace limit of reflection amplitude Ri 2 defined by eq (I4.20p . 
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